AD-A076  204  MASSACHUSETTS  INST  OF  TECH  CAMBRIDGE  GAS  TURBINE  AND— ETC  F/G  21/5 

BOUNOARY  LAYER  AND  WAKE  MODIFICATIONS  TO  COMPRESSOR  DESIGN  SYST— ETC(U 
MAR  79  A  K  SEHRA  F33615-76-C-2118 

UNCLASSIFIED  GT/PDL-144  AFAPL-TR-79-2010  NL 


1 


I  of3 
AO 

407820  ft. 


DDC  .file,  copy  BA0?6204 


AFAPL-TR-79-2010 


.BOUNDARY  LAYER  AND  WAKE  MODIFICATIONS 
TO  COMPRESSOR  DESIGN  SYSTEMS: 

The  Effect  of  Blade-To-Blade  Flow  Variations 
On  The  Mean  Flow  Field  of  a  Transonic  Rotor 


ARUN  K.  SEHRA 


GAS  TURBINE  AND  PLASMA  DYNAMICS  LABORATORY 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
CAMBRIDGE;  MASSACHUSETTS  02139 


MARCH  1979 


FINAL  REPORT 


01  SEPTEMBER  1976  -  31  AUGUST  1978 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  AERO  PROPULSION  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE;  OHIO  45433 

79  11*05  00£ 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  pur¬ 
pose  other  than  in  connection  with  a  definitely  related  Government  procurement 
operation,  the  United  States  Government  thereby  incurs  no  responsibility  nor 
any  obligation  whatsoever;  and  the  fact  that  the  government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications,  or  other 
data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner 
licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any 
rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that 
may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be 
available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


CJLd^t  £-<- 

C.  HERBERT  LAW 
Project  Engineer 


Acting  Chief,  Technology  Branch 


FOR  THE  COMMANDER 


H.  I.  BUSH 


Deputy  Director 
Turbine  Engine  Division 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization  please  notify 
AFAPL/TBX  »  W-PAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a  specific 
document. 

AIR  FORCE/56780/28  SaptamtMr  1979  -  200 


SECURITY^TLASSlFICATION  OF  THIS  PAGE  fl»7i»n  Dmlm  Entered) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

3  .RECIPIENT'S  CATALOG  NUMBER 


1  >  REPORT  DOCUMENTATION  PAGE 


12  GOVT  ACCESSION  NO 


T.Ti.r  ..in  i  i  i.f  - - — — - 

BOUNDARY  LAYER  AND  WAKE  MODIFICATIONS  TO 
COMPRESSOR  DESIGN  SYSTEMS:  The  Effect  of.  Blade' 
To-Blade  Flow  Variations  On  The  Mean  Flow  Field 
of  a  Transonic  Rotor, 


FINAL  REPORT  , 

!}l  Sept— 1976—  31  Aug  1478 

6  PER'FC"HWTTI'S'0W5^R,EPORT-Nm«B*M*— ' " 


GT&PDL  Report  No.  144 


A RUN  K./SEHRA 


10  PROGRAM  ELEMENT  PROJECT.  TASK 


9-  PERFORMING  ORGANIZATION  NAME  AND  AODRESS 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
GAS  TURBINE  &  PLASMA  DYNAMICS  LABORATORY 
CAMBRIDGE,  MA  02139 


AREA  4  WORK  UNIT  NUMBERS 


23074s (A  32 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


AIR  FORCE  AERO  PROPULSION  LABORATORY (TBX) 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433 


March  L979 


14  MONITORING  AGENCY  NAME  »  AOPRESSfll  dllterent  tram  Conlrallin*  Oltice) 


'5.  SECURITY  CLASS,  (of  this  report) 

UNCLASS  I F I  ED 


DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number ) 

TRANSONIC  COMPRESSOR 
BOUNDARY  LAYER 
FLOW  VARIATIONS 
WAKES 


10.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

The  effect  of  b 1 ade-to-b I ade  flow  variations  on  the  mean  flow  field  of  a  hiqhty 
loaded  transonic  axial  flow  compressor  was  investigated.  The  theoretical  approach 
centered  around  modeling  of  three  important  phenomena  associated  with  blade-to- 
blade  flow  fluctuations,  which  control  the  mean  momentum  and  energy  transfer 
processes.  Apparent  stresses  were  introduced  into  the  mean  flow  momentum  equa¬ 
tions  by  pitchwise  averaging.  Loss  concept  of  mean  relative  total  pressure,  due 
to  conversion  of  mean  flow  kinetic  energy  to  the  energy  of  fluctuations,  was 
introduced.  Based  on  this  concept,  mean  rothalpy  and  the  production  of  apparent. 


FORM 
I  JAN  73 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  ffWien  n«r»  Entetn </' 


.UNCLASSlf  IIP. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE(Wf>»n  Data  Enlarad; 


20.  (cont'd) 

entropy  were  defined.  An  expression  for  mean  rotha I py  variation  along  the  stream¬ 
line  was  derived  by  pitchwise  averaging  of  the  energy  eguation.  Mean  flow  egua- 
tions  suitable  for  the  streamline  curvature  computational  scheme  were  developed 
which  include  these  three  effects.  Apparent  stresses,  mean  rothalpy  and  apparent 
entropy  variations  were  calculated  from  measured  velocity  f I uctuat  ions. The  revised! 
streamline  curvature  procedure  was  then  used  to  predict  the  axisymmetric  (peri¬ 
pheral  mean)  flow.  The  agreement  of  rotor  outlet  tangential  velocity  is  excellent 
reproducing  an  unusual  peak  which  is  not  explainable  by  usual  technigues. 

The  results  showed  that  3-0  inviscid  effects  do  not  significantly  modify  the  mean| 
flow,  and  that  the  direct  effect  of  apparent  stresses  on  the  mean  flow  is  also 


sma  I  I  .  jr 
are  due 
to  flow 


Major  deviations  of  the  actual  flow  field  from  the  inviscid  flow  field 
to  the  variations  of  mean  rothalpy  and  production  of  apparent  entropy  duel 
fluctuations,  both  effects  which  are  introduced  here  for  the  first  time. 

\ 


UNCLASS  I F I  ED 


SECURITY  CLASSIFICATION  of  This  PAGE rWiin  Dmtm  Enffd) 


FOREWORD 


This  Final  Scientific  Report  was  prepared  by  the  Gas  Turbine 
and  Plasma  Dynamics  Laboratory  of  the  Department  of  Aeronau¬ 
tics  and  Astronautics,  Massachusetts  Institute  of  Technology, 
under  contract  F33615-76-C-2118  to  the  United  States  Air  Force 
Aero  Propulsion  Laboratory  of  the  Aeronautical  Systems  Division, 
Air  Force  Systems  Command.  C.  Herbert  Law  was  the  AFAPL  Project 
Engineer.  Professor  Jack  L.  Kerrebrock  directed  the  work  with 
the  assistance  of  Arun  K.  Sehra. 

This  Final  Report  describes  the  entire  study  accomplished 
during  the  period  01  September  1976  -  31  August  1978. 


TABLE  OF  CONTENTS 


■ 


Section  Page 

I  INTRODUCTION  1 

1.1  General  Statement  of  the  Problem  1 

1.2  Objectives  of  the  Present  Investigation  3 

1 . 3  Phenomena  Treated  4 

1.4  Method  of  Investigation  4 

II  LITERATURE  SURVEY  11 

2.1  Methods  of  Predicting  the  Flow  Field  12 

2.1.1  Analytical  Methods  12 

2.1.2  Computational  Methods  16 

2.2  Blade  Boundary  Layer  and  Wakes  19 


III  MEAN  FLOW  EQUATIONS  OF  MOTION,  ROTHALPY  VARIATION,  25 


AND  ENTROPY  PRODUCTION 

3.1  Apparent  Stresses  25 

3.2  Mean  Flow  Equations  of  Motion  26 

3.2.1  Non-Bladed  Region  26 

3.2.2  Bladed  Region  29 

3.3  Mean  Flow  Equation  for  Rothalpy  Variation  32 

Along  the  Streamline 

3.4  Mean  Flow  Equation  for  Entropy  Production  40 

3.4.1  Entropy  Production  Associated  with  43 

Flow  Fluctuations 

3.4.2  Entropy  Production  in  the  Rotor  45 

3.4.3  Entropy  Production  due  to  Blade  47 

Passage  Shock 

IV  MEAN  FLOW  EQUATIONS  FOR  STREAMLINE  CURVATURE  49 

COMPUTATIONAL  METHOD 

4.1  Basic  Equation  of  Motion  in  1-Direction  50 

4.2  Equation  of  Motion  for  Computer  Program  54 

4.2.1  Blade  Forces  54 

4.2.2  Equation  for  Meridional  Velocity  58 

Gradient 

4.2.3  Analysis  Type  Momentum  Equation  for  60 

Computer  Program 


4.2.3  Analysis  Type  Momentum  Equation  for  60 
Computer  Program 


1 


TABLE  OF  CONTENTS  (cont'd) 

Section  Page 

V  CRITICAL  ANALYSIS  OF  THE  EXPERIMENTAL  DATA  65 

5.1  Mean  Velocities  67 

5.2  Apparent  Stresses  69 

5.2.1  Apparent  Stresses,  z/c  =0.1  69 

5.2.2  Apparent  Stresses,  z/c  =1.0  71 

5.3  Static  Pressure  71 

5.4  Entropy  Distribution  72 

VI  COMPUTER  PROGRAM  73 

6.1  General  Description  of  Computing  Routine  73 

for  the  Analysis  Program 

6.2  Modified  Streamline  Curvature  Program  79 

6.3  Computation  of  Mean  Flow  Quantities  and  the  82 

Apparent  Stresses  from  the  Experimental  Data 

VII  RESULTS  OF  NUMERICAL  CALCULATIONS  AND  DISCUSSION  85 

7.1  Computational  Mesh  and  Basic  Input  Data  86 

7.2  Mean  Flow  Field  88 

7.2.1  Computed  Flow  with  Cascade  Losses  88 

7.2.2  Computed  Flow  with  Cascade  Losses  89 

and  Apparent  Stresses 

7.2.3  Computed  Flow  with  Cascade  Losses  90 

and  Apparent  Entropy  Production 

7.2.4  Computed  Flow  with  Cascade  Losses  91 

and  Rothalpy  Variation  Along  the 
Streamline 

7.2.5  Computed  Flow  with  Cascade  Losses,  91 

Rothalpy  Vairation,  and  Apparent 

Entropy  Production  Along  the  Streamline 

7.2.6  Computed  Flow  with  Cascade  Losses,  92 

Rothalpy  Variation,  Apparent  Entropy 
Production,  and  Apparent  Stresses 

7.2.7  Computed  Flow  with  Experimental  Losses  93 

VIII  CONCLUSIONS  95 

REFERENCES  99 


vi 


TABLE  OF  CONTENTS  (cont'd) 

Appendix  Page 

A  DERIVATION  OF  RADIAL,  TANGENTIAL  AND  AXIAL  103 

EQUATIONS  OF  MOTION 

B  PITCHWISE  AVERAGING  113 

C  MEAN  FLOW  ENERGY  EQUATION  119 

D  DERIVATION  OF  MEAN  FLOW  EQUATIONS  FOR  MOTION  125 

FOR  STREAMLINE  CURVATURE  TECHNIQUE 


vii 


....  ,•  V  ,f*.  H  J 


LIST  OF  ILLUSTRATIONS 

Figure 

Page 

1 

Rotor  section  cut  a  z  =  constant  plane. 

139 

2 

Passage  shock-wave  configuration. 

140 

3 

Meridional  plane  picture  of  a  passage. 

141 

4 

Blade  forces. 

142 

5 

Experimental  distribution  of  mean  pitch- 
wise  velocity,  0.1  chord  downstream  of 
the  rotor. 

143 

6 

Experimental  Vz  and  V  ,  0.1  chord  down¬ 
stream  of  the  rotor.  r 

144 

7 

Experimental  distribution  of  mean  pitch- 
wise  velocity,  1.0  chord  downstream  of 
the  rotor . 

145 

8 

Experimental  Vz  and  Vr,  1.0  chord  down¬ 
stream  of  the  rotor. 

146 

9 

Apparent  shear  stress  V'V^.,  0.1  chord 
downstream  of  the  rotor. 

147 

10 

Apparent  shear  stress  VZV^ ,  0.1  chord 
downstream  of  the  rotor. 

148 

11 

Apparent  shear  stress  V^.V '  ,  0.1  chord 
downstream  of  the  rotor. 

149 

12 

V'2,  0.1  chord  downstream  of  the  rotor. 

150 

12a 

8vr2  ve 

-Qr~ ■  /  —  ,  0.1  chord  downstream  of  the  rotor. 

151 

13 

2 

Apparent  normal  stress  V'  ,  0.1  chord  down¬ 
stream  of  the  rotor. 

152 

14 

2 

Apparent  normal  stress  V'  ,  0.1  chord  down¬ 
stream  of  the  rotor.  z 

153 

15 

Apparent  shear  stress  V'V' ,  1.0  chord  down¬ 
stream  of  the  rotor.  r  z 

154 

viii 


LIST  OF  ILLUSTRATIONS  (cont'd) 


Figure 

Page 

16 

Apparent  shear  stress  V^,Vg ,  1.0  chord 
downstream  of  the  rotor? 

155 

17 

Apparent  shear  stress  ,  1.0  chord 

downstream  of  the  rotorT 

156 

18 

2 

Apparent  normal  stress  V'  ,  1.0  chord 
downstream  of  the  rotor. r 

157 

19 

2 

Apparent  normal  stress  Vg  ,  1.0  chord 
downstream  of  the  rotor. 

158 

20 

2 

Apparent  normal  stress  ,  1.0  chord 

downstream  of  the  rotor. 

159 

21 

P  /P.  ,  0.1  chord  downstream  of  the  rotor. 

S2  V 

160 

22 

Mean  experimental  entropy,  0.1  chord  down¬ 
stream  of  the  rotor. 

161 

23 

Mean  experimental  entropy,  1.0  chord  down¬ 
stream  of  the  rotor. 

162 

24 

Overall  logic  of  the  original  streamline- 
curvature  program. 

163 

25 

Interaction  of  subroutines  CORDTA ,  SPLINE, 
ROEN1 ,  ROEN2,  and  WATRAN  with  the  main 
program. 

164 

26 

Adjusted  experimental  entropy,  0.1  chord 
downstream  of  the  rotor. 

165 

27 

Layout  of  the  computational  mesh. 

166 

28a 

Basic  input  data  for  computation. 

167 

28b 

Apparent  stress  data. 

168 

29 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

169 

30 

Mean  axial  and  radial  velocities,  0.1  chord 
downstream  of  the  rotor. 

170 

ix 


MW 


LIST  OF  ILLUSTRATIONS  (cont’d) 


Figure 

Page 

31 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

171 

32 

Mean  axial  and  radial  velocities,  1.0  chord 
downstream  of  the  rotor. 

172 

33 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

173 

34 

Mean  axial  and  radial  velocities,  0.1  chord 
downstream  of  the  rotor. 

174 

35 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

175 

36 

Mean  axial  and  radial  velocities,  1.0  chord 
downstream  of  the  rotor. 

176 

37 

_2 

Apparent  stresses  (VQ/r) ,  0.1  chord  down¬ 
stream  of  the  rotor. 

177 

38 

g -F  gg 

Computed  M  and  T(^)app,  station  4. 

178 

39 

Calculated  distribution  of  9I/9m  and 

T  3S/9m,  station  5. 

179 

40 

Apparent  entropy  rise  across  the  rotor. 

180 

41 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

181 

42 

Mean  axial  and  radial  velocities,  0.1  chord 
downstream  of  the  rotor. 

182 

43 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

183 

44 

Mean  axial  and  radial  velocities,  1.0  chord 
downstream  of  the  rotor. 

184 

45 

Calculated  mean  rothalpy  change  across  the 
rotor. 

185 

46 

Mean  pitchwise  velocity,  0.1  chord  downstream 
of  the  rotor 

186 

x 


LIST  OF  ILLUSTRATIONS (cont'd) 


Figure 

Page 

47 

Mean  axial  and  radial  velocities,  0.1  chord 
downstream  of  the  rotor. 

187 

48 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

188 

49 

Mean  axial  velocity,  1.0  chord  downstream 
of  the  rotor. 

189 

50 

Mean  radial  velocity,  1.0  chord  downstream 
of  the  rotor. 

190 

51 

Equivalent  entropy  rise  across  the  rotor. 

191 

52 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

192 

53 

Mean  axial  and  radial  velocities,  0.1  chord 
downstream  of  the  rotor. 

193 

54 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

194 

55 

Mean  axial  velocity,  1.0  chord  downstream 
of  the  rotor. 

195 

56 

Mean  radial  velocity,  1.0  chord  downstream 
of  the  rotor. 

196 

57 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

197 

58 

Mean  axial  and  radial  velocities,  0.1  chord 
downstream  of  the  rotor. 

198 

59 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

199 

60 

Mean  axial  velocity,  1.0  chord  downstream 
of  the  rotor. 

200 

61 

Mean  radial  velocity,  1.0  chord  downstream 
of  the  rotor. 

201 

xi 


LIST  OF  ILLUSTRATIONS  (cont'd) 


Figure 

Page 

62 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

202 

63 

Mean  pitchwise  velocity,  0.1  chord  down¬ 
stream  of  the  rotor. 

203 

64 

Mean  axial  velocity,  0.1  chord  downstream 
of  the  rotor. 

204 

65 

Mean  radial  velocity,  0.1  chord  downstream 
of  the  rotor. 

205 

66 

Mean  pitchwise  velocity,  1.0  chord  down¬ 
stream  of  the  rotor. 

206 

67 

Mean  axial  velocity,  1.0  chord  downstream 
of  the  rotor. 

207 

68 

Mean  radial  velocity,  1.0  chord  downstream 
of  the  rotor. 

208 

xii 


LIST  OF  SYMBOLS 


radial,  pitchwise  and  axial  directions 

v- 

absolute  velocity 

w 

relative  velocity 

time  average  of  velocity  with  respect  to  the 
absolute  frame  of  reference,  or  pitchwise 
average  with  respect  to  the  frame  of  reference 
fixed  to  the  blades 

,<w> 

ensemble  average  with  respect  to  the  absolute 
frame  of .reference,  or  time  average  with  re¬ 
spect  to  the  frame  of  reference  fixed  to  the 
blades 

/  /  .  / 
Vtw,  f> 

fluctuating  part  of  velocity,  pressure 

f  /? 

density 

h  'h 

static  and  total  pressures 

h  - 

reference  pressure 

r,r6 

static  temperature 

Tt 

total  temperature 

• t 

time 

Mean  relative  flow  angle  defined  by 

I'm 

Mean  blade  surface  angle  defined  by  the  mean 
value  of  (  )  for  the  suction  and  pres¬ 

sure  surfaces  of  a  blade 

£ 

blade  lean  angle  in  1-©  plane 

£  7r  V 

Mean  blade  surface  angle  defined  by  the  mean 
value  of  (  -fc  1  for  the  suction  and  pitch' 

wise  surfaces  of  the  blade 

on. 

meridional  streamline  direction 

l 

direction  of  computation  in  the  meridional 
plane  -  outward  from  the  hub 

xiii 


3C'I£I 


LIST  OF  SYMBOLS  (cont'd) 


J  rothalpy  defined  by  1=  h  +  ^ 

^  static  enthalpy 

1*1  total  enthalpy 

_Q  angular  velocity  of  the  rotor 

cjr  total  absolute  velocity  defined  by 

=-A C+W+Vk' 

p  blade  force 

FLtf  Eqv  t-  v  viscous  body  force  components  per  unit  mass 
1  •  in  the  r,0  and  z  directions 

^  radius  of  curvature 

p t  ft  angle  between  the  meridional  and  axial 

directions 

Y  angle  between  the  direction  of  computation  (1) 

and  the  radial  direction,  or  ratio  of  specific 
heats 

'71  direction  normal  to  m-direction  in  the  meridi¬ 

onal  plane 

entropy 

velocity  of  sound 

X  T0v  apparent  stress  terms  appearing  in  the  radial, 
'<p«j  '  pitchwise _and  axial  momentum  equations  __ 

Subscripts  W*-  32  1  *  %  ** 

radial,  pitchwise  and  axial  components 

/yrv  meridional  component,  mean  value 

£  component  in  the  direction  of  computation 

*1  component  perpendicular  to  m-direction 

• 

L  inviscid  component 

IP-  viscous  component 

S  pertaining  to  suction  surface 

b  pertaining  to  pressure  surface 


xiv 


SECTION  I 


INTRODUCTION 

1 . 1  General  Statement  of  the  Problem 

One  of  the  basic  problems  in  achieving  efficient  designs 
of  high  work  compressor  stages  is  the  lack  of  an  adequate 
technique  for  predicting  the  flow  field  and  losses  through  a 
turbomachine  stage.  Current  design  systems  are  based  on  the 
assumption  that  the  flow  can  be  adequately  modeled  by  a  combi¬ 
nation  of  axisymmetric  through  flow  and  two-dimensional 
cascade  f low. ^ ^ ^  whereas  this  procedure  has  proved  quite 
satisfactory  for  lightly  loaded  compressor  stages,  it  is  found 

to  give  a  considerably  different  flow  field  from  that  observed 

2 

xn  highly  loaded  stages.  The  deviation  is  attributed  to 
three-dimensional  effects  which  become  significant  in  a  highly 
loaded  system.  Some  of  the  important  phenomena  responsible 
for  three-dimensional  effects  are  blade  passage  shocks, 
the  inviscid  pressure  field  of  the  blades  resulting  in  the 
distortion  of  stream  surface,  span-wise  flows  in  the  blade 
boundary  layer  and  wakes,  secondary  flows  due  to  the  end  wall 
boundary  layer,  and  shock  boundary  layer  interaction.  The  first 
two  are  inviscid  effects  and  their  proper  treatment  would  require 
a  full  three-dimensional  transonic  inviscid  computation. 

The  last  three  effects  are  intrinsically  viscous  in  origin,  as 
they  involve  boundary  layer  growth  on  the  blade  surfaces  and 
the  end  walls.  The  evolution  of  wakes  and  their  subsequent 
interaction  with  the  adjacent  blade  row  makes  the  problem  of 
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predicting  the  flow  field  more  complex.  Further  complexities 

are  introduced  by  the  three-dimensional  nature  of  the  rotor  blade 
boundary  layer  coupled  with  shock  boundary  layer  interaction,  and 
secondary  flow  arising  from  the  end  wall  boundary  layer  and  tip 
leakage.  An  accurate  estimate  of  viscous  effects  would  require 
a  detailed  computation  of  the  three-dimensional  blade  surface 
boundary  layer  and  end  wall  flow. 

Current  analytical  and  computational  methods  for  prediction 

of  the  flow  field  are  not  advanced  enough  to  account  for  all 

the  coupled  three-dimensional  effects  present  in  a  highly  loaded 

system.  At  best  one  can  account  for  each  individual  effect  by 

isolating  it  from  the  other  three-dimensional  effects.  Although 

2  3 

recent  advances  in  3-D  inviscid  and  multi  2-D  computational 
schemes  are  encouraging,  their  use  is  restricted  to  lightly 
loaded  systems  since  no  capability  exists  to  account  for  the 
above  mentioned  viscous  effects.  It  will  be  shown  later  in  this 
report  that  3-D  inviscid  effects  alone  do  not  produce  significant 
changes  in  the  mean  flow  field  of  the  compressor,  rather  3-D 
viscous  effects  are  mainly  responsible  for  deviation  of  the  mean 
flow  field  from  that  predicted  by  the  usual  axisymmetric  theories. 

One  important  step  in  prediction  of  the  flow  field  of  a 
compressor  is  inclusion  of  the  effect  of  blade  to  blade  (or  pitch- 
wise)  flow  variations  on  mean  (tangentially  averaged)  flow  field 
of  the  compressor.  This  is  the  goal  of  this  investigation. 
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1 . 2  Objectives  of  the  Present  investigation 


More  specifically,  the  present  investigation  has  the 
following  objectives: 

a)  To  determine  the  effect  of  pitchwise  flow  variations 
on  the  mean  flow  field  of  an  axial  compressor,  in 
particular,  to. determine  how  pitchwise  flow  varia¬ 
tions  effect  the  circumferential  mean  momentum  and 
energy  transport  processes. 

b)  To  derive  the  equations  of  mean  flow  for  the  stream¬ 
line  curvature  computational  scheme  implemented  in 
the  computer  program  of  Reference  (1) ,  including 
the  effects  of  blade  to  blade  flow  variations. 

c)  To  carry  out  a  critical  analysis  of  the  time  resol¬ 
ved  experimental  data  of  the  MIT  Blowdown  Compressor 
the  purpose  being  to  assess  the  magnitudes  of  vari¬ 
ous  three-dimensional  effects  (discussed  in  Section 
1.1  of  this  chapter)  in  a  highly  loaded  transonic 
compressor. 

d)  To  modify  the  computer  program  of  Reference  (1)  in 
order  to  account  for  the  effects  of  pitchwise  flow 
variations  on  the  mean  flow  field  and  then  to  com¬ 
pute  the  mean  flow  field  using  the  velocity  correla¬ 
tions  as  estimated  from  the  experimental  data. 


1.3  Phenomena  Treated 

In  order  to  introduce  the  nonaxisynunetric  effects  into 
the  axisymmetric  computation,  the  following  phenomena  are 
treated: 

a)  Apparent  stresses  introduced  in  the  mean  flow  field 
due  to  pitchwise  averaging  of  the  equations  of 
motion. 

b)  Streamwise  variation  of  mean  Rothalpy  due  to  energy 
transfer  across  the  stream  surfaces  by  apparent 
stresses,  and  by  conversion  of  mean  flow  kinetic 
energy  to  energy  of  fluctuations. 

c)  Apparent  entropy  production  due  to  dissipation 
associated  with  apparent  stresses,  and  due  to  the 
production  of  fluctuation  energy. 

1.4  Method  of  investigation 

It  is  well  known  that  the  blade  surface  boundary  layer 
and  wakes  are  primarily  responsible  for  the  pitchwise  flow 
variations  within  and  downstream  of  a  rotor  blade  row.  one 
way  of  accounting  for  the  pitchwise  or  blade  to  blade  flow 
variations  is  to  treat  them  as  flow  fluctuations  around  the 
circumferentially  averaged  flow.  With  respect  to  an  abso¬ 
lute  frame  of  reference  these  will  be  temporal  fluctuations 
around  the  mean  which  is  a  time  mean.  These  fluctuations 
can  therefore  be  treated  as  analogous  to  turbulent  flow 
fluctuations  and  a  quasi-steady  approximation  can  be  applied 
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to  the  flow  field  downstream  of  a  blade  row.  However,  if 
the  frame  of  reference  is  fixed  to  the  blades,  as  is  the 
normal  practice  for  flow  analysis  within  the  blade  row, 
these  blade  to  blade  flow  variations  represent  spatial  fluc¬ 
tuations  around  the  mean  which  is  a  spatial  mean  or  mean  in 
rO.  By  analogy  to  turbulent  stresses  (Reynolds  stresses) 
these  large  scale  fluctuations  will  also  introduce  apparent 
stresses,  and  the  momentum  and  energy  transfer  associated 
with  these  stresses  will  modify  the  mean  flow  field  of  the 
compressor. 

Following  the  above  mentioned  concept,  mean  flow  equa¬ 
tions  of  motion  for  a  frame  of  reference  fixed  to  the  blades 
are  first  derived  (Section  3.2).  This  is  achieved  by  carry¬ 
ing  out  a  pitchwise  average  from  the  pressure  surface  of  one 
blade  to  the  suction  surface  of  the  adjacent  blade.  An 
integral  representation  for  averaging  any  quantity  is  given 
by  Os 

J-J—  / 

£s-©/>  J 

0j> 

where  0S  and  ©p  represent  the  angular  positions  of  the  suc¬ 
tion  and  pressure  surfaces  of  the  two  adjacent  blades.  The 
averaging  procedure  which  uses  Liebnitz  rule  for  finding  the 
pitchwise  average  of  derivatives  has  been  discussed  by  Smith4 
and  Fujii^.  Pitchwise  averaging  of  the  equations  of  motion 
introduces  two  types  of  additional  terms  into  the  mean  flow 
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equations  of  motion.  One  type  of  term  represents  the  in- 
viscid  blade  forces  while  the  other  type  contains  velocity 
correlations  which  can  be  considered  as  apparent  stresses 
in  the  mean  flow,  it  may  be  mentioned  here  that  existing 
axisymmetric  through  flow  theories  do  include  inviscid  blade 
force  terms  but  assume  that  velocity  correlation  terms  are 
small.  A  detailed  derivation  of  the  pitchwise  averaged 
equations  is  given  in  Appendix  a. 

For  the  purpose  of  flow  field  computation,  the  pitch- 
wise  averaged  equations  of  motion  are  derived  for  the  Stream¬ 
line  curvature  Computational  method.  Normal  practice  for  the 
Streamline  Curvature  scheme  is  to  derive  the  radial  equilib¬ 
rium  equation  (  40  )  in  which  the  axial  derivative  terms  are 
replaced  by  the  meridional  and  radial  derivatives,  and  the 
pressure  gradient  is  replaced  by  the  enthalpy  and  entropy 
gradients  as  in  Crocco's  equation.  However,  for  a  more 
general  case  where  a  computing  station  need  not  be  along  a 
radial  direction,  the  method  suggested  by  wennerstron/  is 
more  suitable.  Streamline  Curvature  equations  based  on  the 
Wennerstrom  procedure  have  the  advantage  over  the  usual  type 
that  a  computing  station  can  be  placed  along  the  blade  lead¬ 
ing  or  trailing  edge.  The  equation  of  motion  in  apparent 
stress  format  and  along  an  arbitrary  direction  in  the  merid¬ 
ional  plane,  is  suitably  modified  for  application  to  the 
computer  program  of  Reference  (1) . 
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The  general  problem  of  prediction  of  the  loss  distribu¬ 
tion  is  important  because  in  the  Crocco  form  of  momentum 
equation,  the  distributions  of  entropy  and  rothalpy 
represent  the  pressure  gradient,  and  therefore  these  have  an 
important  effect  on  the  flow  field.  The  losses  can  be  viewed 
as  consisting  of  a  viscous  part  and  an  inviscid  part.  Here 
viscous  losses  refer  tothe  blade  surface  friction  losses  and 
can  be  approximated  from  the  two-dimensional  cascade  data. 
Because  of  the  gradual  variation  of  viscous  losses  in  the 
spanwise  direction,  these  losses  primarily  effect  the  com¬ 
pressor  efficiency  but  have  little  influence  on  the  spanwise 
flow  distribution.  The  concept  of  inviscid  losses  is  not  as 
straightforward  as  the  viscous  losses,  and  in  fact  is  intro¬ 
duced  to  the  turbomachinery  context  here  for  the  first  time. 
It  refers  to  the  conversion  of  mean  flow  kinetic  energy  to 
the  energy  of  fluctuations,  in  the  case  of  turbulent  flow, 
such  a  conversion  is  normally  known  as  "loss  to  fluctuations" 
and  it  implies  that  the  energy  of  fluctuations  is  an  unavail¬ 
able  form  of  energy  as  far  as  the  mean  flow  is  concerned. 
These  losses  can  be  termed  "apparent  losses"  and  their  effect 
is  felt  through  the  reduction  of  mean  relative  total  pressure 
along  the  streamline.  For  the  particular  case  of  the  Crocco 
form  of  mean  flow  momentum  equation,  these  losses  will  alter 
the  mean  rothalpy  and  mean  entropy  gradient  terms.  It  may  be 
mentioned  here  that  the  energy  of  fluctuations  is  eventually 
converted  to  thermal  energy  by  viscous  dissipation  at  the 
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molecular  level,  resulting  in  the  increase  of  real  entropy. 
Here  it  is  assumed  that  the  presence  of  fluctuations  creates 
an  apparent  increase  in  entropy  equivalent  to  that  resulting 
from  the  viscous  dissipation  of  the  energy  of  fluctuations  to 
heat. 

Modeling  of  inviscid  losses  (Section  III)  is  achieved  by 
defining  a  mean  rothalpy  X  and  an  apparent  entropy  produc¬ 
tion  along  the  streamline. 

As  will  be  discussed  in  Section  III,  the  pitchwise 
averaged  rothalpy  is  defined  as 


(1) 


which  would  result  from  the  averaging  procedure  discussed 
earlier. 

The  production  of  apparent  entropy  along  a  streamline 
is  defined  as 
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where  ^  represents  dissipation  associated  with  apparent 

stresses,  and  /z')  represents  the  production  of  flue¬ 

nt 

tuation  energy  along  the  streamline.  The  concepts  used  in 
arriving  at  the  above  definition  are  discussed  in  Section  III. 


An  equation  for  the  rate  of  change  of  mean  rothalpy 
along  the  streamline  is  derived  by  carrying  out  pitchwise 
averaging  of  the  energy  equation  using  equation  (1), 


The  estimation  of  spanwise  distribution  of  mean  rothalpy 
(or  apparent  entropy)  at  any  axial  station  is  achieved  by 
numerically  integrating  the  equations  for  and 

along  the  meridional  streamline,  the  constant  of  integration 
being  evaluated  from  the  known  value  of  mean  rothalpy  (or 
apparent  entropy)  at  the  upstream  station.  This  completes 
the  theoretical  formulation  for  "apparent  stress  correction" 
to  the  mean  flow  model  based  on  the  combination  of  axisym- 
metric  through  flow  and  two-dimensional  cascade  flow.  A 
detailed  derivation  of  the  equations  for  the  mean  rothalpy 
and  apparent  entropy  gradients  along  the  meridional  stream¬ 
line  is  given  in  Section  III. 


Experimental  data  from  the  MIT  Blowdown  Compressor  is 
used  to  estimate  distributions  of  mean  flow  properties  such 
as  velocities,  apparent  stresses  (velocity  correlations) , 
and  mean  entropy.  All  the  flow  parameters  were  reduced  to 
correspond  to  standard  conditions  at  the  inlet  with  air  as 
the  fluid,  velocity  correlations  are  calculated  for  use  in 
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the  theoretical  formulation  for  the  prediction  of  flow  field 
using  Streamline  Curvature  Computational  technique. 

Finally,  the  Streamline  curvature  program  of  Reference 
(1)  is  modified  to  account  for  the  effect  of  blade  to  blade 
flow  variations.  A  splinefit  approximation  is  adopted  for 
interpolating  the  experimentally  estimated  mean  quantities, 
inputs  to  the  modified  streamline  curvature  program  are  the 
hub  and  casing  geometry,  relative  flow  angle  distribution  at 
the  blade  exit,  viscous  loss  distribution  from  the  two- 
dimensional  cascade  data,  mass  flow,  and  apparent  stress 
distribution. 
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SECTION  II 


LITERATURE  SURVEY 

Past  work  being  reviewed  in  this  section  has  been  divided 
into  two  parts,  one  dealing  with  the  methods  of  predicting  the 
flow  field  of  an  axial  compressor,  and  the  other  dealing  with 
blade  boundary  layers  and  wakes.  The  latter  are  responsible  for 
the  large  pitchwise  flow  variations,  spanwise  flow  and  large 
spanwise  gradients  of  losses. 

The  problem  of  predicting  the  flow  field  through  a  turbo¬ 
machine  stage  can  be  approached  either  as  a  design  (indirect) 
problem  or  as  an  analysis  (direct)  problem.  In  the  design  problem, 
the  distribution  of  angular  momentum  change  or  total  enthalpy 

change)  across  the  blade  row  is  assumed  to  be  known.  The  blade 
profiles  have  to  be  chosen  to  give  the  desired  angular  momentum 
change,  with  the  aid  of  experimental  cascade  data.  The  analysis 
problem  deals  with  prediction  of  the  flow  field  from  the  known 
blade  exit  flow  angles  and  the  losses  through  the  blade  row. 

This  problem  is  considerably  more  complex  than  the  design  problem, 
but  is  more  versatile  since  it  can  be  applied  to  all  operating 
conditions  of  the  turbomachine.  The  success  of  analytical  methods 
has  been  mainly  restricted  to  solving  design  problems.  For  solving 
an  analysis  problem  one  has  normally  to  resort  to  a  computational 
method. 

The  blade  surface  boundary  layer  and  wakes  seem  to  have  received 
very  little  attention  from  turbomachine  analysts.  A  few  semi- 
impirical  efforts  to  predict  the  nature  of  blade  boundary  layers 
and  wakes  have  not  proven  good  enough  for  any  realistic  application. 
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2.1  Methods  of  Predicting  the  Flow  Field 


2.1.1  Analytical  Methods 

During  the  past  four  decades  considerable  effort  has  gone  into  the 

methods  of  flow  analysis.  Early  efforts  were  concentrated 
on  the  flow  analysis  for  simplified  flows.  Starting  with 

p 

pitchline  analysis  ,  the  flow  analysis  techniques  were  soon 

g 

extended  to  simple  radial  equilibrium  .  in  the  radial  equil¬ 
ibrium  theory  one  looks  only  at  the  flow  upstream  and  down¬ 
stream  of  the  blade  row  where  the  radial  velocity  is  assumed 
to  be  zero.  The  theory  also  assumes  a  completely  axially 
symmetric  flow.  The  drawback  of  this  theory  is  that  it  gives 
no  account  of  how  the  changes  in  velocity  distribution  take 
place  across  a  blade  row.  A  major  breakthrough  in  the  flow 
analysis  resulted  from  the  concept  of  'Actuator  disc'  with 
the  aid  of  which  not  only  the  velocity  distribution  across 
a  blade  row  could  be  predicted  but  mutual  interaction  be¬ 
tween  the  blade  rows  could  also  be  accounted  for.  in  the 
actuator  disc  theory,  it  is  imagined  that  the  blade  row, 
moving  or  stationary,  is  replaced  by  a  disc  of  infinitely 
small  axial  thickness,  across  which  a  sudden  change  in 
tangential  velocity  and  vorticity  takes  place.  Radial 
equilibrium  exists  far  upstream  and  far  downstream  of  the 
actuator  disc,  but  not  necessarily  in  between  these  stations. 
The  parallel  assumption  of  axial  symmetry  is  also  made,  as  in 
the  radial  equilibrium  theory.  Analytical  studies  of  the 
actuator  disc  model  have  been  made  by  several  turbomachine 
designers^"®'  ^  but  perhaps  the  most  complete  account 
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employing  the  concept  of  actuator  disc  was  presented  by 

Marble*3.  Although  this  theory  has  been  classified  as  three- 
dimensional  since  it  allows  for  velocities  in  the  radial, 
azimuthal  and  axial  directions,  it  really  permits  flow 
variations  in  two  directions  only,  axial  and  radial. 

In  this  sense,  it  is  really  a  two-dimensional  theory, 
in  this  pseudo  three-dimensional  theory  which  has  come  to  be 
known  as  "axisymmetric  through  flow  theory",  the  effect  of 
blades  is  represented  by  distributed  forces  having  the  same 
integrated  effect  as  the  blades  in  the  radial,  azimuthal  and 
axial  directions.  Marble's  success  in  getting  analytical 
solutions  for  several  types  of  circulation  distributions 
(indirect  problem)  resulted  from  the  fact  that  he  was  able 
to  linearize  the  problem  by  considering  the  three-dimensional 
aspect  of  the  flow  as  a  perturbation  on  the  upstream  flow. 
However,  the  assumption  of  small  perturbation  on  the  upstream 
flow  restricted  the  application  of  this  theory  to  lightly 
loaded  systems  only.  For  highly  loaded  systems  it  was  not 
possible  to  linearize  the  problem  and  therefore  computational 
techniques  were  required  for  analyzing  the  flow  field. 

The  above  mentioned  analytical  methods  for  the  predic¬ 
tion  of  turbomachine  flow  have  no  provision  for  estimating 
losses  and  flow  deviation  angles.  This  necessitated  a 
parallel  effort  for  obtaining  the  loss  distribution  and 
deviation  angles  with  the  aid  of  experimental  studies  since 
the  theoretical  methods  failed  to  give  realistic  results. 
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These  experimental  studies  could  be  generalized  only  for 
two-dimensional  cascade  flows  in  which  the  characteristics 
of  through  flow  are  represented  by  the  entering  and  leaving 
flow  angles,  by  a  mean  stream  surface  slope  along  which  the 
cascade  lies,  and  by  the  change  of  axial  mass  flow  density 
between  the  inlet  and  outlet.  Thus  a  combination  of  axis- 
ymmetric  through  flow  and  two-dimensional  cascade  flow  be¬ 
came  an  attractive  technique  for  predicting  the  mean  flow 
field  of  an  axial  compressor.  This  concept  was  generalized 
by  Marble^  for  lightly  loaded  systems  and  was  considered  a 
major  breakthrough  in  solving  a  design  problem.  However, 
with  the  introduction  of  highly  loaded  compressor  stages  for 
modern  high  thrust /weight  ratio  aircraft  engines,  this  line¬ 
arized  theory  was  not  considered  very  attractive  for  designing 
an  axial  compressor. 

An  inviscid  3-D  flow  model  was  achieved  by  McCune  and 
Hawthorne  in  1976.  In  the  case  of  axisymmetric  through  flow 
one  visualizes  the  trailing  vorticity  arising  from  the  nonuni¬ 
formity  of  blade  loading  to  be  smeared  uniformly  in  the  pitch- 
wise  direction.  On  the  other  hand,  if  it  is  assumed  that  the 
trailing  vorticity  is  concentrated  in  discrete  sheets,  it  is 
possible  to  evaluate  the  pitchwise  flow  variation  from  the 

induced  velocity  field  of  the  trailing  vortex  sheets.  McCune 
15 

and  Hawthorne  succeeded  in  using  this  concept  in  their 
three-dimensional  inviscid,  incompressible  theory  for  a 
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rectilinear  cascade.  They  replaced  the  blade  row  with  a  set 

of  lifting  lines  of  bound  vorticity,  and  assumed  that  sheets 

of  vorticity  from  the  trailing  edge  of  each  blade  are  con- 

vected  approximately  by  the  mean  flow  (Beltrami  Flow) .  in 

order  to  apply  this  theory  to  a  highly  loaded  blade  row, 

mean  flow  equations  of  motion  were  not  linearized  and  the 

solution  for  mean  flow  was  obtained  by  using  an  iterative 

schema.  Several  subsequent  improvements  were  introduced  in 

this  theory  to  make  it  more  versatile  in  its  application. 

Morton16  applied  this  theory  to  compressible  flow  through  a 

highly  loaded,  rectilinear  cascade.  Cheng1^  applied  this 

theory  to  incompressible  flow  through  a  highly  loaded  annular 

cascade.  All  the  above  cases  were  restricted  to  nearly  free 

vortex  flow  for  which  the  spanwise  variation  in  loading  is 

18 

not  large.  This  restriction  was  removed  by  Adebayo  and 

the  theory  was  made  applicable  to  a  highly  loaded  blade  row 

with  arbitrary  spanwise  load  distribution.  The  theory  was 

made  possible  by  the  concept  of  treating  the  three-dimensional 

aspect  of  the  mean  flow  field  as  perturbation  about  the  axis- 

ymmetric  mean  flow  instead  of  about  the  upstream  flow,  as  is 

the  case  with  other  linearized  theories.  Therefore,  this 

theory  though  still  ’'linearized''  is  linearized  about  the 
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fully  nonlinear  mean  flow.  Tan  extended  McCune's  three- 
dimensional  theory  to  vorticity  arising  from  viscous  wakes. 
These  three-dimensional  theories  for  predicting  pitchwise 
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flow  variations  downstream  of  a  rotor  are  a  major  step  for¬ 
ward  in  the  analytical  methods.  However,  since  the  knowledge 
of  mean  Ve  distribution  is  required  a  priori,  these  theo¬ 
ries  are  restricted  to  'design  problem'. 

2.1.2  Computational  Methods 

With  the  advent  of  high  speed  digital  computers,  the 

possibility  of  predicting  the  mean  flow  field  of  a  highly 

20 

loaded  compressor  stage  appeared  bright.  Wu  was  probably 
first  to  formulate  equations  for  computing  the  complete 
three-dimensional  flow,  in  his  general  theory  the  equations 
of  fluid  flow  are  satisfied  on  intersecting  families  of 
streamsurfaces,  the  complete  three-dimensional  flow  being 
obtained  by  an  iterative  process  between  the  solutions  for 
the  flow  on  the  two  sets  of  surfaces.  The  two  sets  of 
streamsurfaces  are  blade  to  blade  surfaces  and  s2  surfaces 
which  pass  through  the  blade  row.  However,  Wu's  general 
theory  did  not  find  much  application  and  axisymmetric  coun¬ 
terpart  of  this  theory  (Wu's  through  flow  theory)  was  consid¬ 
ered  more  suitable  for  computing  the  flow  field.  Wu's  through 
flow  theory  is  similar  to  the  general  theory,  but  the  equa¬ 
tions  of  flow  are  solved  for  the  mean  S2  streamsurfaces  on 
which  the  flow  and  fluid  state  are  regarded  as  average  values 
for  the  flow  within  the  blade  passage.  For  a  multistage 
machine,  the  time  dependence  of  the  flow  is  removed  by 
treating  the  through  flow  solution  as  an  axially  symmetric 
flow  for  the  duct  region  between  each  pair  of  blade  rows. 
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During  the  past  decade  two  numerical  methods  based  on 
the  combination  of  axisymmetric  through  flow  and  two-dimen¬ 
sional  cascade  flow  became  the  guidelines  for  most  of  the 

computer  programs.  These  methods  have  come  to  be  known  as 

4  2122 

the  Streamline  curvature  technique^'  '  an(j  the  Matrix 

23 

through  flow  method  .  in  the  streamline  curvature  computa¬ 
tional  procedure,  a  radial  equilibrium  equation  is  solved 
simultaneously  with  mass  conservation  in  the  meridional 
plane.  The  equations  are  written  in  terms  of  the  meridional 
streamline  slope  and  curvature,  which  are  fixed  initially  by 
assuming  uniform  flow  through  the  duct,  and  are  iteratively 
revised  after  solving  the  flow  equations.  The  Matrix  through 
flow  method  solves  the  flow  equations  on  the  s2  surface  of 
Wu’s  through  flow  theory.  Here  the  governing  equation  is 
transformed  into  a  Poisson  equation  for  the  stream  function. 
The  resulting  matrix  equation  is  solved  by  using  a  suitable 
ite  rat i ve  scheme . 

More  recently  at  the  MIT  Gas  Turbine  Laboratory, 

24  2 

Oliver  and  Sparis  ,  and  Thompkins  have  successfully  devel¬ 
oped  and  afplied  a  fully  three-dimensional  inviscid  scheme 
to  a  transonic  compressor  rotor.  Although  this  scheme  is  a 
big  step  in  evaluation  of  three-dimensional  flow  of  a  highly 

loaded  system,  its  utility  is  limited  because  it  does  not 

25 

include  viscous  effects.  Novak  and  Hearsey  have  developed 
a  multi  two-dimensional  technique  for  computing  the  three- 
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dimensional  flow  field.  The  method  is  based  on  the  iterative 
coupling  of  streamline  curvature  solutions  on  Wu's  Si  and  S2 
surfaces  where  is  a  concentric  blade  to  blade  surface  and 
the  S2  surface  represents  a  mean  stream  sheet  defined  by  stream 
function/'s  0.5.  The  main  limitation  of  this  method  is  that 
it  cannot  be  applied  to  a  system  where  large  blade  to  blade 
flow  variations  exist  as  is  the  case  in  the  blade  surface 
boundary  layer. 

Most  of  the  well  known  computational  schemes  use  finite 

difference  methods  for  computing  the  flow  field.  Efforts 

are  being  made  by  a  few  workers  to  apply  finite  element 

computational  methods  to  the  turbomachine  flow  field  but  it 

is  difficult  to  assess  their  effectiveness  at  this  early 

stage.  A  computational  scheme  based  on  the  finite  element 

method  which  examines  the  contributions  of  pitchwise  flow 

variations  on  the  axisymmetric  through  flow,  was  proposed  by 
26 

Hirsch  in  1975  .  The  basic  approach  adopted  by  Hirsch  for 

pitchwise  averaging  of  the  equation  of  motion  was  similar  to 

the  one  undertaken  for  this  investigation,  but  his  apparent 

stress  model  is  unrealistic  because  he  used  the  2-D  isolated 

27 

wake  model  of  Silverstein  et  al. t  for  a  wake  structure 
which  is  highly  three-dimensional,  in  addition,  his  model 
did  not  include  the  concept  of  mean  rothalpy  and  apparent 
entropy,  vhich  will  be  seen  to  be  of  major  impact. 
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It  may  be  mentioned  here  that  all  the  computational 
schemes  discussed  here  depend  on  external  sources  (such  as 
cascade  data)  for  loss  distribution  through  the  blade  row. 

2 . 2  Blade  Boundary  Layer  and  Wakes 

As  mentioned  earlier,  the  viscous  effects  are  mainly 
responsible  for  the  large  blade  to  blade  or  pitchwise  flow 
variation  in  a  highly  loaded  blade  row.  Due  to  the  combined 
action  of  increased  centrifugal  force  in  the  blade  boundary 
layer  and  coriolis  forces  arising  from  the  blade  surface 
curvature,  large  spanwise  flows  are  set  up  in  the  blade 
boundary  layer.  This  makes  more  complex  the  problem  of 
estimating  the  pitchwise  flow  variation.  Further,  since  the 
spanwise  loss  distribution  is  significantly  altered  due  to 
the  spanwise  flows  in  the  blade  boundary  layer,  a  detailed 
computation  of  three-dimensional  boundary  layer  is  needed.  But 
the  existing  knowledge  of  the  computational  methods  suitable 
for  the  present  day  computer  hardware,  are  not  advanced 
enough  to  handle  such  a  problem.  Therefore,  a  virtual  vacuum 
seems  to  prevail  as  far  as  the  theoretical  modeling  of  blade 
boundary  layer  is  concerned. 

The  first  systematic  effort  for  computing  the  three- 

dimensional  boundary  layer  to  a  turbomachine  blade  row  seems 

28  29 

to  have  been  by  Mager  in  1954,  who  adopted  Gruschwitz 

type  of  coordinate  system  for  the  three-dimensional  integral 

momentum  boundary  layer  prediction.  He,  however,  assumed  the 
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blade  surface  to  be  plane  and  introduced  several  other 
simplifications  to  achieve  the  integral  momentum  solution. 
Subsequent  effort  by  several  workers  resulted  in  many 
improvements  and  more  recently  Hules30  was  able  to  formulate 
integral  momentum  equations  for  a  three-dimensional  compres¬ 
sible  boundary  layer  in  a  rotating  variable-pitch  helical 
coordinate  system,  ideal  for  describing  arbitrary  shaped 
turbomachine  blades.  He  computed  the  boundary  layer  charac¬ 
teristic  parameters  for  a  turbine  rotor  without  any  stream- 
wise  pressure  gradient.  Application  of  this  technique  to 
an  axial  compressor  rotor  with  large  adverse  pressure  gradient 
is  still  to  be  tried. 

With  respect  to  the  blade  wakes,  the  major  effort  in  the 

past  has  been  restricted  to  study  of  the  wake  structure  of 

27 

isolated  airfoils.  Silverstein,  Katzoff,  and  Bullevant 
were  the  first  to  investigate  the  near  and  far  wakes  of  an 
isolated  airfoil  in  1939.  They  used  a  theoretical  treatment 
put  forth  by  Prandtl  to  generalize  their  experimental  results 
to  predict  the  important  wake  parameters  in  terms  of  the  dis¬ 
tance  behind  the  trailing  edge  and  the  profile  drag  coeffi- 
.  31 

cient.  Preston  and  Sweeting  in  1945  carried  out  a  system¬ 
atic  investigation  of  the  characteristics  of  the  wake  behind 
an  isolated  airfoil  and  observed  that  a  similarity  in  mean 
velocity  profile  exists  close  behind  the  airfoil  and  that 
the  wake  centerline  velocity  recovered  to  about  80%  of  the 
free  stream  velocity  in  a  quarter  chord  length  from  the 
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in  1953  to 
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trailing  edge.  These  observations  led  Spence 
give  a  general  expression  for  the  wake  centerline  velocity 
as  a  function  of  distance  behind  the  trailing  edge  of  the 
airfoil  and  the  wake  edge  velocity.  No  attempt,  however, 
seems  to  have  been  made  to  arrive  at  a  generalized  theore¬ 
tical  formulation  of  wake  structure  as  a  function  of  the 
physical  characteristics  of  an  airfoil  or  its  loading. 

Very  recently  analytical  and  experimental  investigations  of 
the  characteristics  of  3-D  turbulent  wakes  downstream  of  a 

turbomachine  rotor  have  been  carried  out  by  Raj  and 
32 

Lakshminarayana  .  By  introducing  the  Reynolds  stress  terms 
in  the  equations  of  motion  expressed  in  streamwise  orthogonal 
coordinate  system,  they  tried  to  develop  an  approximate  quasi- 
three-dimemional  turbulent  wake  model  for  a  turbomachine  rotor. 
Their  conclusion  that  the  wake  velocity  defect  almost  com¬ 
pletely  vanishes  very  near  the  blade  trailing  edge,  does 

not  agree  with  the  experimental  data  from  MIT  Blowdown 
2 

Compressor  which  indicates  that  the  velocity  defect  can 
even  grow  with  distance  downstream  of  the  blade  trailing 
edge.  Since  the  rotor  used  by  Raj  and  Lakshminarayana  was 
a  non-lifting  type,  the  validity  of  their  model  for  a  highly 
loaded  compressor  rotor  has  still  to  be  confirmed.  During 
the  past  few  years,  the  MIT  Gas  Turbine  Laboratory  has  been 
conducting  systematic  experimental  and  theoretical  studies 
of  the  behavior  of  wakes  behind  a  highly  loaded  compressor 
rotor.  By  application  of  a  general  small  disturbance  theory. 
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Kerrebrock  has  shown  that  the  wakes  are  not  purely  con- 

vected  in  a  strongly  swirling  flow,  but  can  exhibit  wave 

35 

like  propagating  behavior.  Kerrebrock  and  Mikolajczak 

carried  out  a  kinematic  analysis  of  the  interaction  of  rotor 

wakes  with  stator  and  predicted  the  pitchwise  distribution 

of  stagnation  temperature  and  pressure  at  stator  exit  by 

.  27 

using  the  two-dimensional  wake  model  of  Silverstein  et  al.y 

in  the  generalized  form  as  given  by  Kemp  and  Seers^6.  By 

using  a  fast  response  strain  gauge  probe,  Thompkins  and 
37 

Kerrebrock  successfully  obtained  wake  velocity  distributions 
at  two  axial  locations  downstream  of  a  highly  loaded  transonic 
compressor.  As  mentioned  earlier,  one  of  the  important 
findings  was  the  persistence  of  wake  velocity  defect  down¬ 
stream  of  the  rotor.  Measurements  made  at  a  distance  of  1.0 
chord  behind  the  trailing  edge  indicated  an  increased  velo¬ 
city  defect  at  certain  radial  locations. 

Hence,  it  is  seen  that  most  of  the  previous  investiga¬ 
tions  on  the  flow  modeling  for  axial  compressors  have 
neglected  the  effect  of  pitchwise  flow  variations  while 

2  6 

estimating  the  mean  flow  field.  The  only  attempt  by  Hirsch 
in  trying  to  assess  the  effect  of  pitchwise  flow  variations 
was  inconclusive  since  an  unrealistic  two-dimens xonal  wake 
model  was  adopted  for  modeling  the  velocity  correlations. 

The  author  has  not  come  across  any  publication  which  examines 
the  possibility  of  the  existence  of  inviscid  or  apparent 
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losses.  These  losses  may  have  oreater  effect  on  the  mean 

flow  than  the  viscous  losses. 

It  is  further  seen  that  a  virtual  vacuum  appears  to 

exist  as  far  as  the  study  of  three-dimensional  blade  boundary 

layer  and  wakes  is  concerned.  The  only  detailed  investiga- 

30 

tion  of  three-dimensional  blade  boundary  layer  has  been 

restricted  to  integral  momentum  methods  for  turbomachine 

blades  without  any  pressure  gradient.  Similarly,  the  only 

33 

known  flow  model  for  three-dimensional  wakes  does  not 
appear  to  be  useful  in  its  present  form  for  a  highly  loaded 
compressor  blade  row. 

It  may  be  concluded  from  the  study  of  literature  pre¬ 
sented  in  this  chapter  that  the  existing  knowledge  of  analy¬ 
tical  and  computational  techniques  is  not  advanced  enough 
to  solve  the  problem  of  predicting  the  flow  field  of  a 
highly  loaded  transonic  axial  flow  compressor.  However,  the 
past  work  is  very  useful  in  understanding  several  important 
concepts  and  the  limitations  of  various  analytical  and  com¬ 
putational  techniques.  The  availability  of  the  time  resolved 
experimental  data  downstream  of  the  MIT  Blowdown  Compressor 
rotor  is  an  opportunity  to  assess  the  order  of  magnitude  of 
various  three-dimensional  affects  (discussed  in  Section  i)  . 
This  data  is  of  vital  importance  in  studying  the  effect  of 
pitchwise  flow  variations  on  the  mean  flow  of  an  axial  com¬ 
pressor,  which  is  the  objective  of  this  investigation. 
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SECTION  III 


MEAN  FLOW  EQUATIONS  OF  MOTION,  ROTHALPY  VARIATION  AND  ENTROPY 

PRODUCTION 

3.1  Apparent  Stresses 

It  is  well  known  that  temporal  fluctuations  in  the 

velocity  field  of  a  flow  with  spatial  gradients  influence 

the  mean  motion  in  such  a  way  that  the  mean  flow  exhibits 

am  apparent,  increase  in  resistance  to  deformation,  in  other 

words  these  fluctuations  introduce  "apparent  stresses",  thus 

augumenting  the  viscous  stresses.  For  turbulent  flow  these 

apparent  stresses  are  commonly  known  as  Reynolds  Stresses. 

The  direct  application  of  this  idea  to  the  flow  field  behind 

a  turbomachine  blade  row  would  consist  in  time  averaging  the 

three-dimensional  equations  of  motion  written  in  a  coordinate 

system  fixed  to  the  blade  row  (stationary  for  a  stator  and 

rotating  for  a  rotor) ,  whereupon  Reynolds  stresses  arise  as  a 

result  of  turbulent  fluctuations  (or  temporal  fluctuations) in 

that  coordinate  system.  This  approach  has  been  followed  by 

33 

Raj  and  Lakshminarayana  in  their  attempt  to  study  the 
characteristic  behavior  of  three-dimensional  wakes  downstream 
of  a  rotor. 

Here  the  concept  of  apparent  stresses  is  used  in  a 
different  way.  in  order  to  account  for  the  effect  of  blade 
to  blade  or  pitchwise  flow  variations,  it  is  proposed  to 
consider  them  as  flow  fluctuations  around  the  mean  defined  as 
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axisymmatric  mean.  With  respect  to  the  absolute  frame  of 
reference,  these  pitchwise  fluctuations  would  become  temporal 
fluctuations  around  the  mean  which  is  a  time  mean,  and  analo¬ 
gous  to  the  turbulent  flow,  these  fluctuations  would  introduce 
apparent  stresses  in  the  mean  flow  field.  However,  if  the 
frame  of  reference  is  fixed  to  the  blade  row,  these  pitchwise 
flow  variations  represent  spatial  fluctuations  around  the 
mean  which  is  a  spatial  mean  or  mean  in  r©.  Here  the  apparent 
stresses  would  appear  in  the  equations  of  motion  as  a  result 
of  pitchwise  averaging  and  can  play  a  significant  role  in  the 
transfer  of  mean  momentum  and  energy. 

3.2  Mean  Flow  Equations  of  Motion 

3.2.1  Non-B laded  Region 

Normally,  the  procedure  for  arriving  at  the  governing 
equations  in  Reynolds  stress  format  is  for  a  rotor  exit  flow 
field,  identical  to  the  classical  procedure  described  in 
Reference  (38)  for  unsteady  turbulent  flows. 

It  consists  in: 

a)  Combining  the  continuity  equation  with  the  momentum 
equation. 

b)  Expressing  each  variable  as  the  sum  of  a  mean  plus 
a  fluctuating  part,  for  example  the  velocity. 

U  =  v  +  V 

where  _  0 

—  CA-cm  J— 

7 _ ^  oo  T 
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and  T  is  large  enough  so  that  y  is  independent 

of  T. 

c)  Averaging  the  momentum  equations  in  time  over  the 
interval T,  thus  eliminating  all  terms  linear  in 
the  disturbance  quantities,  while  retaining  those 
which  are  quadratic. 

It  should  be  noted  at  this  point  that  if  the  rotor 
exit  flow  is  fully  periodic  in.  blade  passing,  then  T  need 
be  only  one  blade  passing  period.  But  the  procedure  admits 
the  possibility  of  dealing  with  flows  which  are  unsteady 
in  rotor  coordinates  by  simply  increasing  T  to  average  out 
such  fluctuations. 

If  density  fluctuations  are  neglected,  the  mass  con¬ 
servation  equation  for  the  mean  flow  remains  unchanged  but 
the  equations  of  motion  and  energy  will  have  additional 
terms  containing  correlations  of  fluctuating  quantities. 
However,  if  density  fluctuations  are  not  small,  mean  mass 
conservation  equation  will  also  have  additional  terms 
containing  correlations  of  fluctuating  quantities. 

The  present  analysis  will  be  restricted  to  the  case  of 
small  density  fluctuations.  It  is  further  assumed  that  the 
rotor  exit  flow  is  periodic  in  blade  passing,  which  means 
in  the  absence  of  turbulent  fluctuations  the  temporal 
average  of  unsteady  flow  with  respect  to  the  absolute  frame 
of  reference  is  the  same  as  the  pitchwise  average  with 
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respect  to  the  moving  frame  of  reference.  We  can  thus 
represent  the  effect  of  pitchwise  flow  variations  in  the 
non-bladed  region  by  time  averaging  the  equations  of  motion 
written  in  the  absolute  frame  of  reference.  Mean  mass  and 
momentum  conservation  equations  in  the  absolute  cylindrical 
frame  of  reference  are: 


Mass  Conservation 


_L  2(W)  .  2(%t)  _  Q 


(4) 


Momentum  conservation 
r-direction  2 

iSt  2S  .,14  ^ 


Q-direction 
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z-direction 


Here  the  correlations  of  fluctuating  velocities  repre¬ 
sent  the  effect  of  pitchwise  flow  variations  on  the  axisym.- 
metric  mean  flow  and  should  not  be  confused  with  turbulent 
fluctuations.  Analogous  to  Reynold  Stresses,  these  correla¬ 
tions  can  be  considered  as  apparent  normal  and  shear  stresses. 
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normal  stress  since  it  represents  the  transport  of  radial 
momentum  in  the  radial  direction  due  to  the  above  mentioned 
fluctuations.  Similarly,  a  term  such  as  —  $  '\J\r 

can  be  considered  an  apparent  shear  stress  and  represents  the 
transport  of  pitchwise  momentum  through  a  surface  normal  to  the 
radial  direction,  and  are  the  average  vis¬ 

cous  forces  per  unit  mass  in  the  r ,  0,  and  z  directions. 

3.2.2  B laded  Region 

The  method  of  Reference  (38)  for  deriving  mean  flow 

equations  of  motion  is  not  applicable  within  the  blade  row 

because  of  the  presence  of  blades  across  which  a  jump  in 

flow  properties  takes  place.  The  correct  procedure  for  the 

bladed  region  is  to  write  down  the  equations  of  motion  in 

the  relative  frame  of  reference  and  carry  out  spatial  averaging 

in  the  circumferential  direction  from  the  pressure  surface 

of  one  blade  to  the  suction  surface  of  the  adjacent  blade. 

39 

Ruden  appears  to  have  been  first  to  discuss  the  use  of 

pitchwise  averaged  equations  in  turbomachine  aerodynamics. 

4  40 

The  subject  was  treated  in  greater  detail  by  Smith  ,  Horlock 
and  Fujii5.  This  averaging  procedure  introduces  two  types 
of  additional  terms  in  the  mean  flow  equations  of  motion. 

First,  there  are  those  which  contain  the  effect  of  pressure 
jump  across  the  blade,  while  the  second  are  the  correlations 
of  the  fluctuating  quantities.  The  former  can  be  interpre¬ 
ted  to  represent  inviscid  blade  forces  which  can  be  expressed 
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in  terms  of  the  qradients  of  angular  momentum  by  using  the 
procedure  of  Reference  (7) .  It  may  be  mentioned  here  that 
axisymmetric  through  flow  theory  does  include  inviscid  blade 
force  terms  but  assumes  that  velocity  correlation  terms  are 
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The  pitchwise  or  circumferential  average  of  any  quantity 

i  ,0  2)  is  given  by  the  integral 

9S 

where  0S  and  0p  represent  the  angular  positions  of  the 
suction  and  pressure  surfaces  of  the  adjacent  blades. 

It  is  also  known  from  Liebnitz  Rule  that  integral  of  a 
derivative  is  given  by 

/  ct&  =  iL/*9  +  i  _  /  2& 

4),  **  a*  4  l9) 

Using  equations  (8)  and  (9)  ,  pitchwise  averages  of 

the  derivatives  of  pressure  and  the  products  of  velocities 
can  be  evaluated.  For  example,  pitchwise  averages  of 
and  are  approximately  given  by  the  following 

relations 


2 k  =  +  (Ps  -  6 
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(11) 


+  ^1^2  ■?>  wiwi  3> 
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where  <5~  is  the  blade  lean  angle  and  2 7[)^  _ 

The  term  y^b }  represents  the  radial 

(  ITT?'/ N) 

component  of  inviscid  blade  force.  A  detailed  derivation  of 
equations  (10)  and  (11)  is  given  in  Appendix  b. 

It  is  now  possible  to  obtain  mean  flow  equations  of 
motion  with  the  aid  of  equations  ( io)  and  (11)  .  These  are 
r-direction  momentum  equation 
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9-direction 
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Here  j  Fg-L  and  are  the  r,©,  and  z  components  of 

inviscid  blade  force. 

It  may  be  mentioned  here  that  it  is  possible  to  repre¬ 
sent  the  effect  of  turbulence  explicitly  in  the  above  equa¬ 
tions  by  first  taking  the  time  average  of  the  flow  at  various 
pitchwise  locations  and  then  carrying  out  a  pitchwise  average 
of  the  time  averaged  equations  of  motion.  Appendix  A  gives 
the  details  of  this  procedure. 

Equations  (12)  ,  (13)  ,  and  (14)  can  be  transformed 

to  the  absolute  frame  of  reference  by  putting  ^ 

W©  =  Vg  -  ■fc.Jl 

The  change  in  the  frame  of  reference  does  not  effect  the 
magnitudes  of  the  apparent  stresses.  ;  et<.. 

3.3  Mean  Flow  Equation  for  Rothalpy  Variation  Along  the 
Streamline 

An  important  assumption  in  the  existing  numerical 

methods  is  that  rothalpy  remains  constant  along  the  stream- 

T)  T 

line.  The  Euler  Turbine  equation  .  —  O  ,  which 

X>t 

represents  the  conservation  of  rothalpy  along  the  streamline 
is  applicable  only  when  the  flow  is  steady  and  viscous  ef¬ 
fects  are  not  significant,  it  was  indicated  in  section  1.3 
of  this  thesis  that  the  conversion  of  mean  flow  kinetic 
energy  to  energy  of  fluctuations  will  result  in  the  decrease 
of  mean  relative  total  pressure,  in  order  to  better  under¬ 
stand  this  concept  and  the  correct  meaning  of  the  mean 
relative  total  pressure,  it  is  assumed  that  the  flow  is 


32 


incompressible  so  that  the  relative  total  pressure  can  be 
defined  as 


Pu  - 


b+  £w 

'  2 


(15) 


Pitchwise  averaging  of  this  equation  as  per  the  procedure 
outlined  in  Appendix  B  gives 


fh  =  f>  +  ^  +  SJl" 


(16) 


Here  it  has  been  assumed  that  the  pitchwise  variation  of 
density  is  small  and  the  region  is  outside  the  blade  row. 

Equation  (16)  indicates  that  conversion  of  the  mean 

_ 2 

flow  kinetic  energy  )  to  the  energy  of  fluctuations 

-TV  * 

(fW  )  does  not  change  the  value  of  the  mean  relative  total 
Z 

pressure.  This,  however,  assumes  that  the  energy  of  fluc¬ 
tuations  is  recoverable  when  the  fluid  is  stagnated,  in 
keeping  with  the  physical  postulate  that  the  energy  of  fluc¬ 
tuations  is  not  recoverable  by  the  mean  flow,  we  define  the 
mean  relative  stagnation  pressure  ffe  by 


(17) 
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In  this  example  it  is  assumed  that  the  energy  fed  into 
fluctuations  comes  from  the  mean  flow  kinetic  energy  and  not 
from  the  heat  energy  or  the  energy  related  to  the  pressure 
work. 


Consistent  with  the  concept  of  mean  total  pressure,  the 
mean  total  rothalpy  is  defined  as 


l  =  A  + 

2.  2. 


(18) 


In  order  to  derive  an  expression  for  the  rate  of  change 
of  mean  rothalpy  along  the  mean  streamline,  pitchwise  avera¬ 
ging  of  the  energy  equation  is  carried  out  with  mean  rothalpy 
defined  by  the  above  expression.  The  procedure  involves 
equating  the  mean  rate  of  gain  of  internal  and  kinetic  energy 
to  the  mean  rate  of  work  from  the  pressure  forces,  centri¬ 
fugal  force,  and  the  apparent  stresses.  A  detailed  deriva¬ 
tion  of  ?  I  in  terms  of  the  fluctuating  quantities 

3>t 

is  given  in  Appendix  C.  Equation  (c-22)  from  the  appendix  is 


+ 


(C-22 ) 
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In  deriving  this  equation,  the  work  associated  with 
viscous  stresses  has  been  assumed  to  be  small  as  compared 
to  that  due  to  the  apparent  stresses.  Such  an  assumption 
is  justified  in  the  non-bladed  region,  but  within  the  blade 
row  where  the  viscous  effects  due  to  the  blade  surface 
friction  are  significant,  it  will  be  more  realistic  to  in¬ 
clude  the  work  associated  with  viscous  stresses.  Equation 
(C-22)  for  viscous  flow  will  become 


VI  VII  VIII 


where 


J +  Vc6  w0 

L2^  - - - -j  * 

-f  We  J 


(20) 
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and 


V-f*’ 


2lVa 


For  an  axisymmetric  case 
equation  (19)  will  become 


£  =  H4»  2. 

-2>£  3»v 


(21) 


and  the 


=  .  ^Cfcv/ ^  gj  staple  +  CafZgj-Z  £ 

*  }+=  f-g'+tez  +  R*  ]  ” 

^  J  k$l  xt  5  a  s  gk  j 

+  ==  (  ;  £m/)  l+-Wf& 

i  2>t  3*  J  yy 


(22) 


We  will  now  examine  different  terms  of  equation  (19) 
to  get  a  better  understanding  of  the  energy  transport  pro¬ 
cesses  responsible  for  the  variation  of  mean  rothalpy  along 
the  streamline. 

Term  I  -  j-  j fa  fa  )  +  WeZe  1-WzJgJ  - 

represents  the  work  associated  with  normal  and  shear  apparent 
stresses  as  given  in  equations  (C-15)  and  (C— 16 )  of  Appendix  C. 
These  terms  can  be  interpreted  to  represent  the  work  done  by 
the  forces  due  to  apparent  stresses  working  at  the  mean  velocities. 
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Term  II  =  —  —  ^  '*■  represents  the  rate  of  change 

X>k 

of  the  kinetic  energy  of  fluctuations  along  the  mean  meridi¬ 
onal  streamline.  It  is  seen  that  an  increase  in  kinetic 
energy  of  fluctuations  along  the  mean  streamline  will  result 
in  a  decrease  of  mean  rothalpy,  and  this  is  consistent  with 
the  physical  concept  discussed  in  the  preceding  paragraphs. 

Term  III  =  -  P  t  -  (Wr  .  w'  2t'  ) 


can  be  interpreted  as  work  done  by  the  velocity  fluctuations 
in  correlation  with  pressure  fluctuations,  it  can  in  fact 
be  compared  to  acoustic  energy  transport. 


Terms  IV  and  V  are  interpreted  as  the  production  of  the 
energy  of  fluctuating  pressure  due  to  the  density  gradient. 

Term  vi  can  be  written  as 


and  the  individual  terms  in  the  right  hand  side  of  this 
equation  can  be  expressed  as 


Here  the  first  term  on  the  right  hand  side  is  interpreted  as 
the  spatial  transport  of  the  "kinetic  energy  of  fluctuations" 
due  to  velocity  fluctuations.  The  second  term  represents 
the  production  of  fluctuation  energy  due  to  spatial  variation 
of  fluctuating  velocity. 


Term  VIII  which  represents  the  effect  of  molecular 
viscosity  is  responsible  for  the  spatial  transport  of  kinetic 
energy  due  to  viscous  shear  stresses.  This  spatial  transport 
process  is  always  accompanied  by  viscous  dissipation  resul¬ 
ting  in  loss  of  mean  relative  total  pressure,  it  may  be 
noted  that  in  the  case  of  viscous  dissipation,  the  loss  of 
relative  total  pressure  occurs  due  to  the  conversion  of 
kinetic  and  pressure  energy  to  the  heat  energy.  On  the 
other  hand  the  loss  associated  with  fluctuations  is  due  to 
the  conversion  of  mean  flow  kinetic  energy  to  the  energy  of 
fluctuations. 

Equation  (19)  in  its  entirety  is  not  tractable  pre¬ 
sently  because  Terms  III  through  VIII  cannot  be  estimated 
from  the  available  experimental  data,  in  order  to  estimate 
these  terms  it  is  necessary  to  carry  out  a  local  analysis  of 
blade  to  blade  flow  variations  so  that  all  the  three- 
dimensional  effects  are  properly  accounted  for.  It  may  be 
then  possible  to  arrive  at  a  suitable  flow  model  for  pre¬ 
dicting  blade  to  blade  flow  distribution  required  for 
evaluating  Terms  III  through  VIII.  For  the  purpose  of  the 
present  analysis,  these  terms  will  be  neglected.  That  such 
a  neglect  is  acceptable  can  be  argued  since  these  terms  are 
not  as  large  as  Terms  I  and  II.  For  example.  Term  VIII, 
which  represents  work  associated  with  viscous  shear  stresses, 
will  be  small,  at  least  outside  the  blade  row.  if  the 


spatial  derivatives  of  density  are  smaller  than  the  spatial 
derivatives  of  apparent  stresses.  Terms  IV  and  V  can  be 


neglected.  it  may  be  argued  that  the  spatial  derivatives 
of  fluctuating  quantities  are  smaller  than  the  spatial 
derivatives  of  apparent  stresses,  terms  III  and  VII  can 
therefore  be  neglected.  similarly  since  the  fluctuation  in 


velocities  is  expected  to  be  smaller  than  the  mean  velocity, 
?:(*%)  will  be  small  as  compared  to  Cw  A.') 

3>t  3>t 


in  view  of  the  above  argument,  equation  (22) 
becomes 


(23) 
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In  a  more  compact  form  equation  (23) 


can  be  written  as 


H  =  -  ^ +  tar.fi  I&  +  CjfiZ-z  J 
-Jr  I  _  2 

lLrafr  ^ 

3.4  Mean  Flow  Equations  for  Entropy  Production 


An  established  practice  for  the  inviscid  analysis  of 
turbomachine  flow  field  is  to  write  the  equation  of  motion 
in  Crocco  form  by  using  the  Gibbs  relation 


T  Vs  =  Vk  -  Yt 

f  (25) 

The  final  form  of  Crocco' s  equation  suitable  for  flow 
through  the  rotor  becomes 

~V/X  fVxj/j  =  VI  -  T  Vi  (26) 

where  T  is  rothalpy  defined  as 

It  is  assumed  to  remain  constant  along  a  streamline. 

There  is  an  apparent  inconsistency  in  the  Crocco 
equation  since  it  contains  an  entropy  variation  term  with¬ 
out  including  the  source  of  this  entropy  variation.  One 
way  of  interpreting  such  a  situation  is  to  consider  that 
viscous  effects  are  negligible  at  the  point  where  this 
equation  is  being  applied,  but  upstream  of  this  point 
viscous  effects  do  exist  as  a  result  of  which  an  entropy 
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gradient  is  present  at  the  point  of  application.  Fince 
the  viscous  effects  are  ignored  in  this  ecruation,  the 
entropy  gradient  is  in  a  direction  perpendicular  to  the 
streamline.  Therefore,  the  use  of  this  form  of  erocco's 
equation  can  probably  be  justified  for  application  to  a 
region  outside  the  blade  row  if  the  dissipative  process 
occurs  mainly  within  the  blade  row.  Such  a  situation  is 
quite  unrealistic  for  a  highly  loaded  machine  in  which  one 
of  the  major  sources  of  loss  is  the  dissipation  due  to 
turbulent  mixing  downstream  of  a  blade  row. 

41 

In  order  to  overcome  this  inconsistency,  Horlock 
suggested  that  a  dissipative  body  force  term  could  be  in¬ 
troduced  into  the  equation  of  motion  so  that  the  Crocco 
equation  becomes 

W.  x  (Vx  If  )  —  V-T  —  ^7S  —  (27) 

By  taking  the  dot  product  with  W  ,  we  get 

'BS  =  -  W  •  Z>  (28) 

2>t  — — 

This  is  a  streamwise  ecruation  for  entropy  production  and 
if  a  simplified  loss  model  is  adopted  based  on  the  assumption 
of  dissipative  forces  opposing  the  relative  velocity  vector, 
we  can  express  the  dissipative  force  in  terms  of  a  stream- 

wise  entropy  gradient. 
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(29) 


£>  =  T 

'ds 
42 

Bossman  and  Marsh  suggested  an  improvement  over 
Horlock’s  formulation  by  writing  the  equation  of  motion  in  a 
direction  normal  to  the  streamlines  and  thus  managed  to 
eliminate  the  dissipative  force  term  from  the  equation. 
Wennerstrom7  used  equation  (27)  as  a  starting  point  for 
deriving  the  equation  of  motion  for  a  streamline  curvature 
computational  scheme.  In  this  scheme  also  the  equation  of 
motion  is  written  in  a  direction  almost  normal  to  the  stream¬ 
lines  and  therefore  the  dissipative  force  term  can  be  ne¬ 
glected.  It  may,  however,  be  mentioned  here  that  elimination 
of  dissipative  force  term  by  writing  the  momentum  equation 
in  a  direction  perpendicular  to  the  streamlines,  does  not 
provide  any  special  advantage  because  the  knowledge  of  entro¬ 
py  production  along  the  streamline  is  necessary  to  evaluate 
the  entropy  gradient  term  in  equation  (26) .  This  entropy 
gradient  term  plays  an  important  role  in  changing  the  radial 
distribution  of  velocity. 

Although  computer  programs  (1,42)  based  on  the  above 
model  have  been  successfully  applied  to  turbomachine  flow 
problems,  there  is  a  major  discrepancy  in  this  model  which 
renders  it  unsuitable  for  application  to  a  highly  loaded 
system.  This  formulation  assumes  that  the  experimental 
cascade  data  is  an  adequate  means  for  estimation  of  entropy 
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increase  through  the  rotor,  hence  no  provision  exists  to 
account  for  entropy  rise  associated  with  the  flow  fluctuations. 

In  the  next  section  the  concept  of  "apparent  entropy 
production"  resulting  from  the  flow  fluctuations  is  intro¬ 
duced  . 

3.4.1  Entropy  Production  Associated  with  Flow 
Fluctuations 

In  Section  3.3  an  expression  for  the  rate  of  change  of 
mean  rothalpy  along  a  streamline  was  finalized,  it  was 
emphasized  that  the  spatial  transfer  of  fluctuation  energy 
and  the  streamwise  rate  of  change  of  fluctuation  kinetic 
energy  play  an  important  role  in  changing  mean  rothalpy 
(or  mean  relative  total  pressure)  along  the  mean  meridional 
streamline.  The  concept  of  "apparent  entropy  production" 
due  to  the  flow  fluctuations  will  now  be  introduced. 

The  increased  rate  of  viscous  dissipation  due  to  flow 
fluctuations  is  a  well  known  phenomenon.  The  viscous 
dissipation  from  kinetic  energy  to  thermal  energy,  results 
in  the  decrease  of  relative  total  pressure  and  an  increase 
of  entropy.  Based  on  this  concept  it  should  be  expected 
that  a  decrease  in  the  mean  relative  total  pressure  would 
result  in  an  increase  in  the  mean  entropy,  or  the  conversion 
of  mean  flow  kinetic  energy  to  the  energy  of  fluctuations 
will  cause  the  mean  entropy  to  rise,  it  is  proposed  here 
that  the  production  of  fluctuations  at  any  point  in  the  flow 
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field  is  equivalent  to  the  production  of  apparent  entropy 
the re ,  the  amount  of  apparent  entropy  being  equal  to  the 
kinetic  energy  of  fluctuations  divided  by  the  temperature. 

It  is  also  proposed  that  the  gradient  of  this  apparent  en¬ 
tropy  plays  the  same  role  as  the  real  entropy  gradient  in 
reshaping  the  flow  field.  The  term  "apparent  entropy 
production"  refers  to  the  conversion  of  the  mean  flow  energy 
to  the  energy  of  fluctuations,  the  corresponding  loss  in  the 
mean  relative  total  pressure  can  be  termed  "apparent  losses". 


It  was  argued  in  section  3.3  that  the  mean  relative 
total  pressure  for  incompressible  flow  should  be  defined  by 


equation  (17)  ,  reproduced  here 


ft*.  = 


(30) 


goes  as  . 


Using  the  assumption  that  all  the  energy  of  fluctuations 
— —x  _ a 

(  ^  W_  )  comes  from  the  mean  flow  kinetic  energy  (  f  W  ) t 

2.  2- 
it  is  seen  from  the  above  equation  that  the  rate  of  drop  of 

mean  relative  total  pressure 

3>t 

And  since  the  increase  of  entropy  ~7~  -ZLf  goes  as  __  . 

—  r  f  2>t 

we  postulate  that  7”  goes  as  2>  •  Apparent 

2>t  3>t 

entropy  production  along  a  streamline  also  results  from  the 
apparent  dissipation  which  is  analogous  to  the  viscous 
dissipation  and  refers  to t he  irreversible  work  done  by  the 
apparent  stresses.  Terms  of  the  type  ^ 

represent  apparent  dissipation. 
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Based  on  the  above  argument,  it  is  now  possible  to  write 
an  expression  for  the  apparent  entropy  production  along  a 
streamline,  which  is  as  follows 


r  £f£2#>_  x  .  5 

JDt  T>t 

For  the  axisymmetric  case 

—  =  zL  2. 

-Z>£ 

Equation  (31)  therefore  becomes 


7-tL  _ 


=■$« 


gftv'Zj 


Substituting  for  '  we  get 

T sr*  =  -%[■ **  sp  * *?* 

+  +1K  ) 

+  WgMe  25  l  + 

*  %  y 


3.4.2  Entropy  Production  in  the  Rotor 

Current  practice  for  estimation  of  entropy  increase 
through  the  rotor  is  to  use  the  experimental  data  from  two- 
dimensional  cascade  tests.  But  it  is  well  known  that  the 
entropy  increase  across  a  highly  loaded  rotor  is  larger  than 
that  estimated  from  the  results  of  a  corresponding  cascade 
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unit.  This  behavior  is  more  pronounced  in  regions  near  the 
hub  and  the  casing,  and  at  off  design  operating  conditions. 
Recent  experimental  studies  carried  out  at  MIT  have  confirmed 
that  the  cascade  data  does  not  provide  an  answer  to  the  above 
problem  for  a  highly  loaded  axial  transonic  compressor  rotor. 

Within  the  blade  row,  both  viscous  dissipation  and  the 
apparent  entropy  production  associated  with  blade  to  blade 
flow  variations  are  important.  A  proper  a  priori  estimate 
of  the  change  of  entropy  through  the  rotor  would  require  a 
detailed  study  of  blade  surface  and  end  wall  boundary  layers, 
blade  passage  shocks,  secondary  flows  due  to  end  wall  boundary 
layers,  and  other  three-dimensional  effects  discussed  in 
Section  I  of  this  report.  No  such  detailed  calculation  is 
attempted  here.  Rather,  it  is  proposed  to  estimate  the 
entropy  increase  through  the  rotor  by  assuming  that  the  mean 
flow  viscous  losses  can  be  found  from  the  experimental  cas¬ 
cade  data,  and  the  apparent  entropy  increase  associated  with 
blade  to  blade  flow  variations  with  the  aid  of  equation  (33) 
and  measured  blade  to  blade  flow  variations.  The  entropy 
increase  due  to  the  blade  passage  shock  can  be  found  by  the 
method  given  in  Section  3.4.3. 

Thus  the  total  entropy  increase  through  the  rotor  = 

)  app  +  cascade  +  ^  shock. 
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3.4.3  Entropy  Production  Due  to  Blade  Passage  Shock 


in  order  to  arrive  at  a  realistic  estimate  of  entropy 
production  due  to  blade  passage  shock,  it  is  necessary  to 
know  the  exact  location  and  strength  of  the  shock.  The  prob¬ 
lem  is  made  more  complex  due  to  shock-boundary  layer  inter¬ 
action  resulting  in  the  separation  of  turbulent  boundary 
layer.  For  the  purpose  of  estimating  the  entropy  production 
across  the  shock,  it  is  normally  assumed  that  the  shock 
stands  near  the  entrance  of  the  blade  passage (Figure  2), 
striking  the  suction  surface  at  the  point  b  ,  extending  in 
front  of  the  blade  at  point  A,  and  then  bending  back  similar 
to  a  bow  wave,  it  is  then  assumed  that  the  loss  across  the 
'hock  can  be  approximated  by  the  normal  shock  loss  taken  for 
the  average  of  Mach  numbers  at  points  A  and  b.  The  Mach 
number  at  point  a  is  assumed  to  be  equal  to  the  inlet  relative 
Mach  number.  The  Mach  number  at  B  is  estimated  by  using 

Prandtl-Meyer  expansion  on  the  blade  suction  surface.  This 
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procedure  has  been  suggested  by  Miller  et  al. 

in  the  light  of  the  experimental  studies  carried  out  at 

the  MIT  Gas  Turbine  Laboratory,  the  above  procedure  does  not 

appear  to  be  realistic.  Density  measurements  made  across  the 
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blade  passage  shock  using  a  flow  visualization  technique 
indicated  that  the  flow  Mach  number  at  the  point  3  is  not 
higher  than  that  at  the  point  A.  Thus  the  simple  Prandtl- 
Meyer  expansion  is  not  applicable  for  the  three-dimensional 
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flow  existing  in  the  actual  machine.  It  is  in  fact  more 
realistic  to  assume  that  the  mean  relative  Mach  number  at 
the  blade  inlet  represents  the  uniform  Mach  number  ahead  of 
the  blade  passage  shock.  Entropy  rise  across  the  shock  can 
then  be  estimated  by  the  normal  shock  relation  given  in 


Reference  (45) . 

^5=  S.  ZL 

J  (rw;* 


(34) 
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SECTION  IV 


MEAN  FLOW  EQUATIONS  FOR  STREAMLINE 
CURVATURE  COMPUTATIONAL  MODEL 

In  order  to  study  the  effect  of  blade  to  blade  flow 
variations  on  the  mean  flow  field  of  a  highly  loaded  tran¬ 
sonic  compressor,  a  computational  scheme  commonly  known  as 
the  Streamline  Curvature  technique  was  adopted.  it  is  a 
computational  scheme  based  on  the  assumption  of  axially  sym¬ 
metric  flow  field  and  requires  loss  distribution  as  input 
from  external  sources  such  as  cascade  data.  in  its  most 
general  form  the  method  consists  of  writing  a  single  equa¬ 
tion  of  motion  in  the  meridional  plane  along  the  direction  of 
computation.  For  a  non-bladed  region  the  direction  of  com¬ 
putation  is  normally  taken  as  the  radial  direction,  but  for 
flow  field  computation  at  the  inlet  and  outlet  to  a  blade 
row,  the  direction  of  computation  is  made  to  coincide  with 
the  meridional  projection  of  leading  and  trailing  edges 
respectively.  Since  most  of  the  cascade  data  such  as  flow 
angles,  losses  are  known  at  the  blade  trailing  edge  and 
depend  on  the  conditions  at  the  blade  leading  edge,  the 
above  mentioned  procedure  of  writing  the  equation  along  the 
blade  leading  and  trailing  edges  increases  the  accuracy  of 
the  computational  scheme.  This  single  equation  of  motion  is 
suitably  transformed  by  expressing  radial  and  axial  deriva¬ 
tives  in  terms  of  the  derivatives  in  the  meridional  stream- 
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line  direction  and  the  direction  of  computation.  Streamline 
curvature  enters  into  the  equation  through  t he  terms  con¬ 
taining  radius  of  curvature  and  the  slope  of  meridional 
streamlines.  The  continuity  equation  is  used  along  with  the 
momentum  equation  to  fix  the  constant  of  integration. 
References  (1,3)  have  given  a  detailed  description  of  Stream¬ 
line  Curvature  Technique.  The  aim  of  this  chapter  is  to  re¬ 
derive  the  equations  of  motion  for  the  Computer  program  of 
Reference  (1)  so  that  the  effect  of  blade  to  blade  flow  vari¬ 
ations  can  be  accounted  for. 

4. 1  Basic  Equation  of  Motion  in  1-Direction 

In  order  to  derive  the  mean  flow  equation  of  motion  in 
the  Streamline  Curvature  format,  the  method  adopted  is  to 
first  transform  the  mean  flow  equations  of  motion  in  r,  Q, 
and  z  directions  ( equations  D-l ,  D-2,  and  D-3  of  Appendix  D) 
to  equations  of  motion  in  terms  of  the  components  of  vorticity. 
These  equations  are  then  combined  to  give  a  single  vector 
equation  which  is  further  modified  using  Gibbs  relation  to 
get  a  mean  flow  equation  of  Crocco’s  type. 

Equation  (D-16),  a  detailed  derivation  of  which  is  in¬ 
cluded  in  Appendix  D,  has  the  form 

i«  (v,£)  =vi  -  tvs 

+  ( 'Lz+y-t-  e  Taf.+  k2*i)-  Ti  -  fjf  (d-16) 
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in  this  equation,  the  effect  of  blade  to  blade  flow  varia- 

/A  * 

tions  enters  directly  through  the  terms  y ■+  0 
^■21^  '  T^e  terra  which  is  the  blade  force  vector 

originating  from  the  inviscid  pressure  field  of  the  blade, 
also  enters  the  equation  through  the  pitchwise  averaging 
process  as  a  jump  of  pressure  across  the  blade  from  suction 
to  the  pressure  surface.  Although  the  usual  axisymmetric 
through  flow  formulation  does  include  this  term  as  the  dis¬ 
tributed  forces  representing  the  effect  of  blades,  it  assumes 
that  blade  to  blade  flow  variations  are  small.  The  term  fv 
represents  the  mean  viscous  force  in  the  fluid  and  is  mainly 
due  to  the  blade  surface  friction.  Outside  the  blade  row 
Fi  does  not  exist  and  F i>  is  small. 

The  equation  of  motion  in  the  1-direction  is  derived  by 
taking  the  scalar  product  of  the  vector  equation  (D-16)  with 
a  unit  vector  in  1-direction  and  using  transformation  rela¬ 
tions  for  radial  and  axial  derivatives.  The  final  form  of 
Streamline  Curvature  equation  as  derived  in  Appendix  D 


becomes 


(14  2-Ha  Jcr, Oe+r)+£ Cxft+r) 


-  v\fe.  - 2Wo  Jl Gir 

~f~  Y  -h  Y 

—  ( Fni  "f~  )■£*,(>*+*■)  1 
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Before  commenting  on  this  equation,  the  assumptions 
made  in  deriving  this  equation  can  be  summarized  as: 

a)  Steady  state  flow  with  respect  to  a  coordinate 
system  fixed  to  the  blades. 

b)  Pitchwise  variation  of  density  and  temperature 
are  small. 

It  may  be  noted  that  over  bars  have  not  been  used  for  quanti¬ 
ties  which  are  assumed  to  be  pitchwise  invar ient,  such  as 
f,T,  Y  etc.  The  mean  slope  of  the  streamline  is 
tan  S  = 

Equation  (35)  is  the  1-direction  pitchwise  averaged 

_ t 

momentum  equation  for  U*  .  The  Streamline  Curvature  pro¬ 
cedure  which  begins  by  laying  a  mesh  in  the  meridional  plane 
formed  by  the  initially  assumed  streamlines  and  fixed  axial 
stations,  provides  information  on  ^  and  .  information  on 
other  unknown  quantities  is  obtained  as  follows: 

a)  If  the  mean  rothalpy  distribution  is  specified  at 
an  upstream  station,  it  can  be  estimated  at  the 
desired  station  by  numerically  integrating  equation 

(23)  along  the  meridional  streamline,  when  the 
heat  transfer  and  viscous  effects  are  small. 

b)  Entropy  distribution  at  the  desired  station  is 
estimated  by  using  the  procedure  given  in  Section 
3.4.2. 
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c)  For  design  problem,  1*4  distribution  is  specified  at 

the  blade  trailing  edge.  in  the  interblade  region 
it  can  be  assumed  to  remain  constant  along  the 
streamline.  For  analysis  problem  Wq  can  be  re¬ 
placed  by  the  relative  flow  angle  and  meridional 

velocity  (4,  . 

d)  if  the  value  of  Vjn  is  assumed,  temperature 
can  be  estimated  from  rothalpy. 

e)  Apparent  stresses  are  estimated  from  experimental  data. 

It  is  thus  seen  that  equation  (35)  can  be  integrated 
numerically  along  1-direction  using  an  iterative  technique, 
provided  °  and  the  blade  force  terms  can  be  expressed 

in  terms  of  the  known  quantities.  As  mentioned  earlier,  the 
constant  of  integration  is  provided  by  the  use  of  continuity 
equation.  Further,  since  this  1-direction  equation  of  motion 
could  be  derived  without  using  the  0-direction  equation  of 
motion,  one  may  be  misled  to  the  conclusion  that  r  and  z 
direction  apparent  stresses  are  sufficient  to  get  the  final 
form  of  the  equation  of  motion.  This,  however,  is  not  true 
because  0-direction  apparent  stresses  are  indirectly  intro¬ 
duced  through  the  ©-direction  blade  force  term  F&  which 
is  given  by  pitchwise  averaged  0-direction  momentum  equation 
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Eauation  of  Motion  for  Computer  program 


As  seen  in  the  previous  section,  for  the  application  of 
equation  (35)  ,  it  is  necessary  to  express  the  term 

r~d?n, 

and  blade  force  terms  in  terms  of  the  known  quantities. 

This  equation  may  also  require  further  modifications  in 
order  to  make  it  applicable  to  the  analysis  problem.  Hearsey, 
whose  streamline  curvature  program'*',  will  be  modified  for 
this  investigation,  arrived  at  the  desired  equation  by 
evaluating  _ 'wt  with  the  aid  of  continuity  equation,  and 

'&>rt 

for  the  blade  forces,  he  adopted  the  formulation  given  by 
Wennerstrom^ .  Equations  being  derived  in  this  section  are 
kept  identical  to  those  given  in  Reference  (1) ,  so  that  the 
basic  logic  of  Hearsay's  program  is  not  changed.  Detailed 
derivation  of  only  those  equations  which  are  changed  as  a 
result  of  pitchwise  averaging  is  included. 

4.2.1  Blade  forces 


in  equation  (35)  mean  (or  distributed)  force  terms  in 


the  meridional  direction  and  in  a, direction  perpendicular 
the  meridional  direction  (in  the  meridional  plane)  have  been 
separated  into  the  inviscid  blade  forces  (  F-*» *.t-  Fn.1  )  and 

viscous  body  forces  (  )  .  If  the  viscous  body 

forces  are  considered  to  result  from  the  blade  drag  force 
only,  the  problem  of  estimating  these  forces  becomes  rela¬ 
tively  less  cumbersome. 
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The  basic  method  for  evaluating  the  blade  force  terms 
13 

was  given  by  Marble  ,  which  was  later  generalized  by 
Wennerstrom^  for  any  arbitrary  direction  in  the  meridional 
plane. 

If  the  blade  surface  geometry  is  known,  it  is  possible 
to  fix  the  directions  of  Fi  and  F^  .  which  arises  as 

a  result  of  inviscid  pressure  field  of  the  blade,  can  be 
taken  as  perpendicular  to  the  mean  blade  surface,  and  if  the 
mean  flow  velocity  is  assumed  to  be  approximately  tangential 
to  file  mean  blade  surface,  the  inviscid  blade  force  will  be 
perpendicular  to  the  relative  velocity  vector.  For  evalua¬ 
ting  the  viscous  force  ,  it  is  assumed  that  this  repre¬ 

sents  blade  surface  friction  and  it  is  assumed  to  be  acting 
in  a  direction  opposite  to  the  mean  relative  velocity  vector. 
This  in  fact  is  not  a  realistic  assumption  for  a  highly 
loaded  rotor  in  which  the  boundary  layer  fluid  may  have  a 
significant  radial  component  due  to  the  combined  action  of 
centrifugal  and  coriolis  forces.  However,  in  the  absence 
of  exact  knowledge  of  boundary  layer  flow  and  since  the  vis¬ 
cous  force  term  is  not  very  important  in  the  momentum  equa¬ 
tion,  the  above  assumption  was  retained  for  this  analysis. 

Blade  Force  in  the  Meridional  Direction 

With  the  above  assumption  and  for  small  variation  of 
temperature  and  entropy  in  the  pitchwise  direction,  it  can 
be  shown  that  (for  example,  see  Reference  13) 
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(37) 


T M*  _  A*.iy 


9>n 

Equation  (37)  requires  further  investigation  but  in  view  of 
the  above  discxission,  this  relation  has  been  retained  with¬ 
out  any  alterations. 

From  Figure  4  ,  we  have  Fo  Cos  ft  —  —  W  da*i  ft 

Therefore,  equation  (37)  becomes  (36) 


=  ~r?£ 

We  also  find  from  Figure  4 

~LcL,n ft 

But  from  equation  (36) 

5  -  5, 

Combining  the  above  equations,  we  ge'c 

Cyvi,  =  •  -f-  -  -  YC&  CPF1  ^  t  dj  y  -+• 

Equation  (40)  is  valid  in  the  bladed  region  only. 

For  non-bladed  region,  Fql=0  and  F^  = 
can  be  assumed  to  represent  internal  fluid  friction. 
This  is  given  by 

F**  =  -  T  21  Co/ p 


(38) 


(39) 


(40) 


which 


(41) 
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However,  if  the  friction  is  limited  to  the  blade  surface 
only  then  in  the  non-bladed  region  both  and 

will  vanish. 


Blade  Force  in  a  Direction  Perpendicular  to  the  Meridional 
Direction  ’  . . 

Since  is  assumed  to  be  parallel  to  the  meridional 
streamline  ,  RyLV-  =  ® 

•  r*  —  C 

o «  •  71  ~~  'VL  c 

A  detailed  derivation  of  Rti  is  given  in  Reference  (7  )  . 
Combining  equations  (25)  and  (27)  of  Reference  ( 7 ) ,  we  get 
Rn  =  Mnfi  7"  I ^ 2  7^  Cjf 


Cod  <f  Cod  p  Cod  5 


rh  Lf  ’<&n/3  (0+r)  Ga  e 

Cod  5  t  ..  /  ) 

+  Cod  fi/don.  t  /d£c(f +r)  > 


Here  [p-.y  )  has  been  changed  to  {<pd~T  )  since  Y  represents 
blade  trailing  edge  location  in  Figure  (1)  of  Reference  (7) . 

Equation  (42)  can  be  transformed  to  the  apparent  stress 
format  by  using  equation  (36) ,  thus 

Fn  =  T2lJ 

X  tan(0j.r)  CojE+  CoJjsJCn  £  Jtc  Qf+rJ  l  (43) 
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4.2.2  Equation  for  Meridional  velocity  Gradient 


The  gradient  is  eliminated  from  equation  (35) 

9m 

through  the  use  of  continuity  equation. 

Using  the  transformation  relation 


(D-21 


/ 

and  the  geometric  relationship  — —  -  - - — 

9  “W  St 

the  following  relation  can  be  derived  with  some  algebraic 
manipulation 

_  A&rvrf  DU*  .(l- ta*T>  6  ta.n.r)'}Vk 
9 m  <Co4  r  ?z 

-h  ding}  Co<l  4  2ift:a/n^  (44 

^  Cosr  ^ 

IVz 

The  term  - c  is  now  eliminated  with  the  aid  of  continuity 

equation.  Although,  the  axisymmetric  system  represents  a 
blade  free  space,  it  is  expedient  to  recognize  at  this  point 
the  possibility  of  blockage  due  to  blades.  It  is  considered 
as  a  distributed  blockage  to  avoid  a  conflict  with  the  con¬ 
dition  of  axial  symmetry. 

Denoting  blockage  by  ^  ,  we  can  write  open  area  as 

>=  U-V 

For  the  axisymmetric  model  with  distributed  blockage, 
the  mean  flow  continuity  equation  is  written  as 

2 lOtf>Uk) 


This  equation  is  converted  into  the  following  suitable  form 
for  #  ,  with  the  aid  of  transformation  relation  (D-20) 

-23  =  3:21  +  u*  2h  3^. 

S  'Sm  *t.  A  Cad  T(J.-  *] 

_  ’ta**rtan.4_  7%*  +  ^  Cj*j 
f/-tan/tanrj  9m  Cajr(i-tarx^Unr)^ 

_  icLn  c 

(l-ta.n./£cLnrJ>li  (46) 

It  is  seen  from  equation  (46)  that  we  have  still  to 
eliminate  the  meridional  derivative  of  density  from  the 
first  term  in  the  right  hand  side  of  equation  (46) .  it  may 
be  noted  here  that  although  meridional  variation  of  density 
gradient  is  permitted,  the  pitchwise  variation  of  density 
has  been  assumed  to  be  small. 

The  term  is  obtained  from  -t-  by  using 

-f  ***-  /  3?n 

Second  Law  of  Thermodynamics  and  the  equation  of  state. 

is  obtained  by  deriving  mean  flow  equation  of 

>  9m 

motion  in  the  meridional  direction,  a  detailed  derivation 
of  which  is  given  in  Appendix  D. 

Equation  (D-31)  of  Appendix  D  is 

T  ( 4r  <£»/- 

K.J  4  t  j;  7 

^  (47. 
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/ 

Substituting—  — l  from  equation  (47)  into  equation 

Y  'd'm.  _ 

(46)  ,  and  then  substituting  —  from  (46)  into  (44)  , 

we  get  after  some  algebraic  manipulation  an  expression 
for  ,  given  as 

^iKn  Q  /%*)=  &evn0+ta''nY)  _  %»-  ^<*0 

r&r/l  ( £-  kcLn\fita.-nV) 

-  jX  ^  2J* 

^  ,X  3>n-  ^  'S-m, 

-  ■Mpifv-Qdi'nfzS  7*6  2^^  l 

Vm  l  ^  J 

_ z 

~  -  fal/i* /  -’‘’ZzffrlfijJ  (48) 

Mw.  I 

This  is  the  equation  for  the  meridional  velocity  gradient. 

4.2.3  Analys is-Tyoe  Momentum  Equation  for  Computer 
Program 

To  get  the  final  form  of  mean  flow  equation  of  motion 
in  1-direction  suitable  for  computer  program,  and  the 

blade  force  terms  are  eliminated  from  equation  (35)  by 
using  equations  (48),  (40),  and  (43). 


This  gives  after  some  algebraic  manipulation 
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final  form  of  this  equation  in  the  bladed  region,  suitable  for  the  analysis 
lem  can  now  be  arrived  at.  This  is  given  as 


This  equation  will  be  used  for  computing  the  flow  field 
of  the  MIT  Blowdown  Compressor  rotor.  A  few  important  clar¬ 
ifications  about  this  equation  are  as  follows: 

a)  Equation  (50)  is  applicable  for  the  analysis 
problem  within  the  blade  row  only.  For  the  un- 
bladed  region  where  the  relative  flow  angles 
are  not  specified  and  where  all  blade  related 
information  such  as  blade  forces, ^  f  £  are 
nonexistent,  this  equation  is  suitably  modified 

and  it  takes  the  format  suitable  for  design  problem. 

b)  Although  is  not  specified  for  this  equation, 

it  is  found  from  the  knowledge  of  fi  ,  blade  velo¬ 
city  and  the  assumed  value  of 

c)  Whereas  the  mean  axial  and  radial  velocities  have 
been  expressed  in  terms  of  the  mean  meridional 
velocity,  the  velocity  correlations  (apparent 
stresses)  have  been  retained  as  the  products  of 
the  axial  and  the  radial  velocity  fluctuations. 

This  is  because  the  experimental  data  available 
for  testing  this  formulation  could  provide  velocity 
fluctuations  in  the  radial,  axial  and  the  azimuthal 
directions  only.  The  radial  derivative  of  the 
apparent  stresses  has  also  been  retained  for  the 
same  reason.  The  axial  derivative  of  the  apparent 
stresses  can,  however,  be  transformed  to  the  radial 


^SSSBSSSSBS 
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and  the  meridional  derivatives  by  using  the  relation 


d)  The  mean  rothalpy  gradient  ~%T  is  estimated  with 

'it 

the  aid  of  equation  (23)  . 

e)  The  entropy  gradients  and  in  this  equation 

are  estimated  from  the  cascade  data  and  from  the 
entropy  production  model  as  formulated  in  Section  III. 

f)  It  is  also  observed  that  the  effect  of  radial 
direction  apparent  stress  is  greatest  because  the 
effect  of  pitchwise  and  axial  apparent  stresses  are 
scaled  down  due  to  the  small  multiplication  factor. 
This  is  reasonable  since  the  1-direction  momentum 
equation  approximately  represents  an  r-direction 
momentum  equation. 
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SECTION  V 


CRITICAL  ANALYSIS  OF  THE  EXPERIMENTAL  DATA 

The  experimental  data  adopted  for  the  verification  of 

the  flow  model  formulated  in  the  preceding  chapters,  was 

taken  from  the  MIT  Blowdown  compressor  test  facility.  Time 

resolved  data  on  fluctuating  pressures  had  been  recorded  by 
2 

Thompkins  behind  a  highly  loaded  transonic  rotor  with  the 
aid  of  a  four  diaphragm  probe.  The  data  which  was  digitized 
every  5  Jlsecs  provided  the  instantaneous  values  of  the  three 
components  of  flow  Mach  number  t  /#%  )  and  pressure  at 

different  radial  locations  and  at  two  axial  stations.  These 
axial  stations  were  situated  at  0.1  and  1.0  chords  downstream 
of  the  rotor  blade  trailing  edge.  This  data  was  recorded 
with  respect  to  an  absolute  frame  of  reference. 

For  the  estimation  of  instantaneous  flow  velocities  from 
Mach  numbers,  it  is  necessary  to  know  static  temperature 
which  was  not  measured.  The  static  temperature  was  estimate*' 
by  assuming  the  Euler  Turbine  Equation  is  valid  and  can  be 
expressed  in  the  following  form  for  the  estimation  of  total 
temperature  ratio  across  the  rotor. 


Here  subscripts  1,2  and  T  refer  to  rotor  inlet  and  out¬ 
let  conditions,  and  blade  tip  condition. 
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In  order  to  make  the  results  representative  of  usual 
compressor  tests,  the  inlet  conditions  were  standardized  to 
lti  -  52  o-  0*/?  ,  /Jj  =  7  4>jA‘n2  ,  and  g  =  SJJ 

£  'R.  •  The  choice  of  these  values  at  the  inlet 

result  only  in  multiplication  of  all  the  terms  in  the  momen¬ 
tum  equation  by  a  constant  factor,  the  relative  magnitude  of 
these  terms  remaining  unchanged. 

The  procedure  adopted  for  estimating  the  distribution  of 
velocity  correlations  (apparent  stresses)  Vx  t  ^ 

/  /  f  ~~ i 

14  ,  involved  two  steps.  These  are: 

a)  computation  of  the  time  averaged  values  of 
each  component  of  velocity. 

b)  Computation  of  the  products  of  fluctuating 
quantities  at  each  instant  and  then  time 
averaging  the  products  to  obtain  the  apparent 
stresses . 

Two  points  may  be  noted.  Firstly,  the  correlations  of 
the  absolute  velocity  fluctuations  are  identical  to  those  of 
the  relative  velocity  fluctuations.  secondly,  the  time 
averaged  flow  with  respect  to  the  absolute  frame  of  reference 
is  equivalent  to  the  pitchwise  averaged  flow  with  respect  to 
a  frame  of  reference  fixed  to  the  blades. 

Computations  were  made  of  the  following  quantities,  at 
several  different  radial  locations  and  at  two  axial  stations 
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a) 

b) 

c) 

d) 


Mean  velocity  components  1 / 


Apparent  stresses  y  'Zfc'tfc'  ,  ,  24'*, 


/ 


and 


Static  pressure  ratio  /2r2y/^ 

Mean  entropy  change  across  the  rotor. 


in  computing  the  mean  values,  a  number  of  different 
averaging  intervals  were  used.  Averages  over  one  blade 
passing  period  were  first  constructed,  with  the  averaging 
interval  beginning  at  several  different  points  of  the  blade 
passing  period,  if  the  flow  were  truly  periodic  with  blade 
passing,  that  is  steady  in  rotor  coordinates,  the  mean  values 
should  be  independent  of  the  averaging  interval  so  long  as  it 
covered  one  complete  blade  passing  period.  in  fact  the  aver¬ 
ages  varied  strongly w ith  the  location  of  the  interval,  in¬ 
dicating  that  the  flow  is  quite  unsteady  in  the  rotor  coor¬ 
dinates.  Averages  were  then  taken  over  two  and  three  blade 
passing  periods  starting  from  the  first  data  point. 

Finally,  averages  were  taken  over  five  full  blade  pass¬ 
ing  periods,  the  maximum  available  in  digital  form,  and  these 
values  will  be  interpreted  as  true  means.  The  scatter  from 
these  means  will  be  taken  as  a  measure  of  the  unsteadiness  in 
the  rotor  coordinates. 


5 . 1  Mean  Velocities 

The  mean  velocities  and  t&t  at  the  port  5 

location,  which  is  just  0.1  chord  downstream  of  the  rotor 
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trailing  edge  are  shown  in  Figures  5  and  6.  The  mean 

over  five  blade  passing  periods  is  indicated  by  circular 

points  through  which  the  curve  has  been  drawn.  The  scatter 

from  this  mean,  of  the  values  averaged  over  one,  two  and 

three  blade  passing  periods,  is  also  shown.  In  all  the 

figures  from  Figure  5  through  21  ,  plus  (+) ,  square  (  □  ) , 

triangle  )  and  circle  (o)  represent  means  over  one,  two, 

three  and  five  blade  passing  periods  respectively.  In  some 

of  these  figures  symbol's  Q,  Q  and  *  have  also  been  used. 

These  symbols  refer  to  means  over  one  blade  passing  period 

but  starting  at  different  points  of  the  recorded  data. 

It  is  seen  that  the  scatter  is  small  in  the  inner  part  of 

the  annulus,  very  large  at  %  =  0.87,  decreases  for  some 

distance  towards  the  casing  and  then  becomes  quite  large  near 

the  casing.  The  large  scatter  near  =  0.87  which  was 

2 

first  noted  by  Thompkins  ,  could  be  due  to  a  vortex  shed  from 
the  rotor  blade  near  the  termination  point  of  the  shock  on 
the  suction  surface.  This  view  is  supported  by  the  fact 
that  the  variations  in  v*  are  much  smaller  than  those  in 
t and  . 

Figures  7  and  8  give  the  average  values  of  , 

%  and  'Ufc  at  the  port  6  location,  which  is  one  chord 
downstream  of  the  rotor  trailing  edge.  The  scatter  is  com¬ 
paratively  reduced,  and  the  very  large  V&  near  the  tip  has 
decayed.  There  is  also  a  general  upward  shift  in  in  the 

inner  part  of  the  annulus,  and  a  uniform  upward  shift  in  ^ 
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between  port  5  and  port  6.  The  data  at  the  two  locations 
were  derived  from  two  separate  tests,  and  this  shift  is 
tentatively  attributed  to  a  zero  shift  problem  in  the  pressure 
transducers.  Thus,  only  the  shape  of  the  curves  at  the  two 
stations  should  be  compared.  This  radial  variation  of  both 
V©  and  does  not  change  very  much  with  distance  down¬ 

stream. 

The  mean  radial  velocity  shows  an  erratic  behavior 

at  port  5  and  this  is  attributed  to  a  possible  flow  separa¬ 
tion  near  the  hub,  and  shock-boundary  layer  interaction  in 
the  vicinity  of  radius  ratio  0.87.  No  reasonable  explanation 
can  be  given  for  the  observed  steep  rise  of  1^,  near  the 
casing.  This  is  attributed  to  the  flow  disturbance  generated 
by  the  cavity  which  houses  the  probe.  At  port  6,  however, 
the  erratic  behavior  of  has  considerably  diminished. 

5.2  Apparent  Stresses 

5.2.1  Apparent  Stress,  z/  =  0.1 

The  apparent  shear  stresses,  ,  V^V~q  ,  and 

'\y>lVLQ  are  Shown  in  Figures  9,  10,  and  11  for  the  port  5 
location,  immediately  behind  the  rotor.  Here  again  the  five 
blade  passing  period  means  are  shown  by  the  circular  points, 
and  the  scatter  from  this  mean  is  shown  by  other  point 
symbols  (listed  in  Section  5.1).  Several  points  should 

noted.  First,  the  scatter  is  quite  small  at  r/  ^  **  0.78, 

.  irge  at  r/rt  =  0.87,  and  moderate  near  the  hub  and  the 
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casing.  Secondly,  the  apparent  shear  stresses  —  £  Vj' 

_2 

would  be  of  the  order  of  15  lb^  ft  if  the  compressor  had 
operated  in  ambient  air,  and  hence  are  fairly  small  compared 
to  the  pressure.  On  the  other  hand  there  appear  to  be  rather 
steep  gradients  in  the  stresses,  so  their  effect  may  be  sig¬ 
nificant  nevertheless. 


The  apparent  normal  stresses  and  ^  which 

are  plotted  in  Figures  12,  13,  and  14,  are  much  larger 
than  the  apparent  shear  stresses,  and  also  show  larger 
scatter,  particularly  near  -  0.87.  Taking  mean  values 

over  five  blade  passing  periods  as  representative,  it  is 
observed  that  there  are  very  large  gradients  of  and  1}q 

near  this  radius.  For  example,  the  change  of  from 

nearly  zero  to  30,000  ft2/sec2  between  =0.8  and  0.87 

yields  an  effective  pressure  gradient  of  the  order  of  1000 
tty  /ft  per  unit  radial  distance,  or  in  other  terms  half  the 
static  pressure  per  compressor  radius.  An  alternate^  way  of 

1  Vi* 

assessing  the  order  of  magnitude  of  the  term 

7 

is  to  compare  it  with  the  centrifugal  force  which  approx¬ 

imately  balances  the  radial  pressure  gradient.  A  plot  of 
/  Vq.  VS  is  shown  in  Figure  12a. 

From  this  figure  we  find  that  the  magnitude  of  - *  which 

is  an  important  apparent  stress  term  in  the  1-direction 

momentum  equation,  can  become  four  times  that  of  the  centri- 

.  z 


fugal  force  term 


Therefore,  unless  the  effect  of 
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the  radial  normal  apparent  stress  is  nullified  by  the  other 
apparent  stresses,  it  will  play  a  role  comparable  to  the 
centrifugal  force  at  several  radial  locations. 


5.2.2  Apparent  Stresses  z/c  =  1.0 

As  can  be  seen  from  Figures  15,  16,  and  17,  the 
apparent  shear  stresses  do  not  decrease  very  much  between 
z/c  =  0*1  and  z/c  ~  1*0*  but  there  is  a  large  change  in  the 
radial  distribution  of  stresses.  This  is  consistent  with  the 
tentative  conclusion  of  Reference  (37),  that  the  flow  is 
evolving  to  a  rotor  wake  structure  controlled  by  the  swirl 
eigen  modes  of  the  duct.  The  apparent  normal  stresses  which 
are  plotted  in  Figures  18,  19,  and  20  are  again  seen  to  be 
much  larger  than  the  apparent  shear  stresses. 

5.3  Static  Pressure 

The  distribution  of  mean  static  pressure  divided  by  the 
inlet  stagnation  pressure  A*A ,  is  shown  in  Figure  21 

for  z/c  =0.1  and  z/c  =  1.0.  The  general  trend  is  a  decrease 
in  the  downstream  direction,  which  is  consistent  with  the 
increase  in  •  A  low  pressure  at  ’Az.  =  0.87  directly 

U 

behind  the  rotor  is  probably  associated  with  the  strong 

vortex  which  is  conjectured  at  that  radius,  and  which  is 

.T7I 

also  considered  as  an  explanation  for  the  large  1%,  and  Vq 
at  that  radial  location. 
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5 . 4  Entropy  Distribution 

The  distribution  of  mean  entropy  rise  (  )  was  esti¬ 

mated  by  first  finding  the  instantaneous  entropy  rise  A.5 
with  the  aid  of  the  following  equation,  and  then  carrying  out 
the  time  averaging. 


A  plot  of  - —  Iss  't'/tf.  for  port  5  is  shown  in 

Gr  . — 

Figure  22.  The  negative  values  of  ==—  between  the  radius 

Cy 

ratios  of  0.7  and  0.85  are  confusing  since  the  entropy  change 
being  estimated,  has  been  treated  as  a  thermodynamic  property 
which  must  increase  unless  heat  transfer  takes  place  from  the 
midspan  region  towards  the  hub  and  the  casing. 

In  the  absence  of  any  logical  explanation, this  unrealis¬ 
tic  entropy  distribution  was  attributed  to  the  zero  shift  in 
the  pressure  tranducer.  The  corresponding  plot  for  port  6 
which  is  seen  in  Figure  23  ,  does  not  show  this  anomaly,  so 
that  the  zero  shift  problem  is  mainly  associated  with  measure¬ 
ments  at  port  5.  It  is  also  seen  that  towards  the  casing  the 

value  of  for  port  5  suddenly  becomes  very  large,  even 

Gr 

larger  than  that  at  the  corresponding  location  of  port  6. 

This  again  conflicts  with  die  concept  of  entropy  change  along 
a  streamline  and  indicates  that  either  the  use  of  Euler  Turbine 
equation  is  not  justified  or  the  experimental  data  near  the 
casing  is  not  reliable. 
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SECTION  VI 


COMPUTER  PROGRAM 

As  mentioned  in  Section  I,  one  of  the  objectives  of  this 
investigation  is  to  modify  a  computer  program  based  on  the 
Streamline  Curvature  Computational  technique,  so  that  the  effect 
of  blade- to-blade  flow  variations  can  be  accounted  for  while 
computing  the  mean  flew  field  of  an  axial  flow  compressor. 

The  basic  computer  program  chosen  for  the  purpose  has  been 
developed  by  Hearsey,^  and  is  one  of  the  most  versatile  and 
complete  programs  for  application  to  a  design  or  an  analysis 
problem.  Before  describing  the  modifications  introduced  in  this 
program,  it  is  considered  desirable  to  outline  some  of  its 
important  features  and  the  computational  logic  adopted  for  it. 
This  section  gives  a  general  description  of  the  computing  routine 
for  an  analysis  problem,  the  overall  logic  of  Hearsey's  program, 
and  the  modifications  introduced  in  the  basic  program  to  account 
for  blade-to-blade  flow  variations.  A  brief  description  of  the 
logic  used  for  estimating  the  experimental  flow  from  the  data 
of  Ref.  2  has  also  been  included  towards  the  end  of  this  section. 

6 . 1  General  Description  of  Computing  Routine  for  the  Analysis  Problem 

(Not  Including  the  Effect  of  Pitchwise  Flow  Variations) 

6.1.1  Input  Data 

In  order  to  solve  an  analysis  problem  by  the  Streamline 
Curvature  Computational  technique  the  following  input  data  is 
required  for  a  multistage  axial  flow  compressor. 

a)  Axial  location  and  orientation  of  each  computing 

station:  Generally,  a  number  of  computing  stations 
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are  located  at  strategic  points  in  the  flow  and 
preferably,  but  not  necessarily,  close  to  orthogonal 
to  the  local  meridional  streamline  direction. 

Typically,  several  stations  are  situated  upstream 
of  the  first  blade  row.  Additionally,  stations  may 
be  placed  within  the  blade  row. 

b)  Geometric  details  in  respect  of  hub,  casing  and  the  blading. 

c)  Mass  flow  rate  and  rpm. 

d)  Inlet  conditions  such  as  pressure  and  temperature. 

e)  Distribution  of  losses  through  the  blade  row:  If  the 
losses  are  to  be  specified  on  the  basis  of  cascade 
data,  their  distribution  need  not  be  included  in  the 

input  data.  An  additional  routine  can  be  introduced 

« 

in  the  main  program  which  first  calculates  '  D'  factor 
and  then  losses  from  the  inbuilt  curves.  These  inbuilt 
curves  give  losses  as  a  function  of  '  D*  factor  and 
radial  location  of  the  meridional  streamline. 

f)  The  number  of  streamlines  to  be  used  in  the  computation: 
Generally,  the  location  of  streamlines  which  is  com¬ 
puted  by  the  program,  is  based  on  the  fixed  percentage 
of  the  total  mass  flow  through  each  streamline. 

However,  if  it  is  required  to  concentrate  more  stream¬ 
lines  in  any  particular  region,  it  can  be  achieved  by 
prespecifying  different  percentages  of  the  total  mass 
flow  through  the  individual  streamtubes. 

Hearsey's  program  has  a  provision  for  30  axial  computing 
stations  and  21  streamlines. 
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6.1.2  The  Overall  Solution  Procedure 

The  overall  solution  procedure  corresponding  to  the 
input  data  of  Section  6.1.1,  can  be  described  by  the  following 
sequence  of  steps: 

a)  The  first  step  in  the  computing  process  is  to  form 
a  computational  mesh  by  the  intersection  of  the 
fixed  axial  stations  with  the  meridional  streamlines, 
henceforth  referred  to  as  streamlines.  The  initial 
placement  of  the  streamlines  is  achieved  by  dividing 
the  annulus  area  in  accordance  with  the  number  of 
streamtubes. 

b)  Having  chosen  the  initial  locus  of  each  streamline, 
the  program  establishes  the  radius  of  curvature  (\) 
and  the  slope  (  ^  )  of  each  streamline  at  the  com¬ 
puting  mesh  points  by  using  the  standard  finite 
difference  format 
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c)  Now  the  first  iterative  pass  begins.  During  each 
iterative  pass  the  computation  of  the  flow  field 
is  carried  out  at  all  the  stations,  starting  at 
the  first  computing  station  and  proceeding  to  the 
last.  The  position  of  streamlines  during  a  pass 

is  not  changed.  At  the  end  of  each  pass  the  stream¬ 
lines  are  relocated  by  iterative  process,  and  the 
new  position  of  the  streamlines  is  used  for  the 
next  pass.  The  details  of  calculations  involved 
during  any  computing  pass  are  described  in  the  next 
few  steps.  As  an  example,  the  computing  station 
will  be  taken  at  the  exit  of  the  rotor. 

d)  Before  computing  the  flow  field  at  the  rotor  exit 
during  any  pass,  the  flow  field  at  the  upstream 
station  is  computed.  This  means  that  the  values  of 

velocities,  pressure,  entropy,  etc.  are  already 
known  along  all  the  streamlines  at  the  rotor  inlet 
which  can  be  designated  as  (I-l)th  station. 

To  start  with,  distribution  of  Un  is  assumed  at  the 
Ith  station.  These  values  are  normally  taken  as 
those  obtained  from  the  previous  pass  ( (?l-i)tk  pass)  . 
With  the  assumed  values  of  VL  and  the  given  values 
of  jS  at  each  mesh  point,  the  values  of  can 

be  estimated  using  the  relation 

Vq  =  tcvnp  +  'c/l  (56) 
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Now,  from  the  assumption  of  constant  rothalpy 
along  the  streamline,  the  total  enthalpy  at  the 
rotor  exit  can  be  calculated. 

M (i-i)  -  -til  Ve(i-i)  =  M(l)~  fcJl  Vgfr)  <57, 

- — y - - - ' 

known  at  the  rotor  inlet  from  equation  (56) 

e)  With  total  enthalpy  and  total  velocity  known, 
static  enthalpy  and  thus  the  flow  Mach  number  can 
be  immediately  found.  The  entropy  distribution  is 
found  directly  from  the  specified  losses,  or  these 
losses  are  calculated  by  estimating  the  'D'  factor 
and  then  using  one  of  the  "loss  vs.  '  D*  factor" 
curves  built  into  the  computer  program.  Thus  all 
the  quantities  required  for  integrating  the  momentum 
equation  (50)  are  known. 

f)  Using  the  midstreamline  velocity  value  as  a  starting 
point  (i.e.  as  an  integration  constant) ,  the  momen¬ 
tum  equation  (50)  is  solved  in  a  stepwise  fashion, 
giving  new  values  of  velocities  for  the  inner  and 
the  outer  halves  of  the  annulus.  Once  the  new  values 
of  velocities  are  known  at  all  the  mash  points  of 
the  Ith  station,  it  is  possible  to  estimate  the  new 
values  of  static  enthalpy  and  therefore  temperature. 
Finally,  the  density  distribution  is  estimated  by 
estimating  the  pressure  with  the  aid  of  specified 
losses. 
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g)  With  known  density  distribution,  total  mass  flow 
through  the  annulus  is  computed  and  compared  with 
the  given  value  of  mass  flow.  (it  may  be  noted 
that  the  mass  flow  through  the  individual  stream- 
tubes  is  not  compared  at  this  point.)  in  general 
the  computed  mass  flow  will  not  equal  the  actual 
mass  flow.  The  midstreamline  velocity  is  increased 
or  decreased  as  required,  and  the  process  is  repea¬ 
ted  until  continuity  is  satisfied.  It  is  thus  seen 
that  the  simultaneous  solution  of  momentum  and 
continuity  equations  is  an  iterative  process. 

h)  Continuity  is  now  applied  for  each  streamtube. 

In  general  the  streamtubes  will  not  contain  their 

proper  percentages  of  mass  flow.  The  computer 
determines,  and  stores  an  error  function  for  stream¬ 
line  location. 

i)  The  information  regarding  velocities,  angles,  Mach 
number,  etc., at  the  rotor  exit  is  now  available. 

The  flow  field  at  the  next  downstream  station  can 
therefore  be  computed.  In  the  interblade  region, 
angular  momentum  (or  total  enthalpy)  and  entropy 

are  assumed  to  remain  constant  along  the  streamlines. 

j)  The  above  process  is  continued  to  the  last  computing 
station.  The  previously  determined  streamline 
location  error  functions  are  now  simultaneously 
applied  to  redefine  the  entire  streamline  pattern 
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at  all  the  stations,  and  the  next  iterative  pass 
begins.  Convergence  is  achieved  when  the  change 
in  the  flow  field  is  within  a  specified  limit  for 
two  consecutive  passes. 

The  logic  of  Hearsey's  Computer  Program  is  given  in 
Figure  24.  Here  it  has  been  assumed  that  the  losses  have 
been  given  as  input,  a  detailed  description  of  the  program 
is  given  in  References  (1,3,6). 

6 . 2  Modified  Streamline  Curvature  Program 

in  order  to  modify  the  original  streamline  curvature 
program  to  account  for  blade  to  blade  flow  variations,  five 
additional  subroutines  are  introduced.  The  basic  logic  of 
the  original  program  was  retained  without  any  alteration. 

The  process  of  modification  is  centered  around  estimating 
Rothalpy  variation  and  Entropy  production  along  the  meridional 
streamline  due  to  pitchwise  flow  fluctuations,  and  intro¬ 
ducing  apparent  stress  terms  in  the  1-direction  equation  of 
motion.  Necessary  input  variables  are  introduced  to  study 
the  effect  of  the  above  phenomena  individually  or  in  any 
desired  combination.  These  subroutines  can  only  be  called 
from  "overlay  UD0308"  corresponding  to  NEQN=0  (Reference  1) . 

No  provision  is  made  for  calling  these  subroutines  from 

"overlay  UD0326".  Terms  representing  blade  thickness  (i^,^? 

Ix-i  3*,- 

etc.)  appearing  in  equation  (50)  were  also  not  introduced 
in  these  subroutines.  Necessary  comment  cards  giving  the 
basic  function  of  the  subroutines,  list  of  variables,  and 


79 


other  relevant  details  have  been  introduced.  some  of  the 


details  in  respect  of  these  subroutines  are  as  follows: 

a)  Subroutine  CQRDTA 

This  subroutine  reads  the  velocity  correlation  data 
at  all  the  stations  and  finds  the  cubic  splinefit 
coefficients  with  the  aid  of  subroutine  SPLINE. 

These  coefficients  are  stored  for  estimation  of 
velocity  correlations  and  their  derivatives  at  the 
mesh  points.  Velocity  correlation  data  which  these 
subroutines  read  is  generated  in  a  separate  program 
EXPTL  FORTRAN  from  the  time  resolved  experimental 
data.  This  subroutine  is  called  by  specifying  input 
variable  NCOR  =  1. 

b)  Subroutine  SPLINE 

This  subroutine  calculates  the  coefficients  for 
cubic  splinefit  from  the  input  data  on  velocity 
correlations. 

c)  Subroutine  ROEN  1 

This  subroutine  estimates  the  distributions  of 
Rothalpy  and  Entropy  associated  with  blade  to 
blade  flow  fluctuations.  This  is  done  by  inte- 

Sratln9  <  H  >app  and  <T  |£  >app  the 

meridional  streamline,  the  constant  of  integration 

being  evaluated  from  the  known  value  of  rothalpy 

(or  apparent  Entropy)  at  the  upstream  station. 
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Values  of  ,  T ,  and  7"i£  are  also  esti- 

di  3^  30 

mated  for  use  in  the  momentum  eauation  (50) . 


This  subroutine  is  called  by  specifying  input 
variable  NROEN  =  1. 
d)  Subroutine  roen 

This  subroutine  calculates  (  2L  )  _  and  (7~^£  ) 

2^  aPP  -dTK  aPP 

with  the  aid  of  equations  (23)  and  (33)  respectively. 
Apparent  stress  terms  and  their  radial  derivatives 
are  computed  by  using  the  cubic  splinefit  coefficients 
generated  in  subroutine  CORDTA. 

Axial  derivatives  of  apparent  stresses  are  calculated  by  first 
estimating  the  meridional  derivatives  and  then  using  the  re¬ 
lation  / 


Meridional  derivative  of  any  quantity  can  be  calculated 
by  using  the  standard  finite  difference  approximation 


5Z  )  2.0 

'  ^1+1  -  **2  ~  J /  (59) 

If  input  variable  NROTH  =  1,  (  )  is  calculated, 

r 

otherwise  (  ___  ) app  is  taken  equal  to  zero.  Similarly,  if 

input  variable  NENTR  is  not  equal  to  1,  (  7^1  )  is  taken 

*2?*,  aPP 

to  be  zero. 
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e)  Subroutine  WATRAN 


This  subroutine  calculates  apparent  stress  terms 
appearing  in  the  pitchwise  averaged  1-direction 
momentum  equation  (50) .  Cubic  splinefit  coefficients 
generated  in  subroutine  CORDTA  are  used  to  estimate 
apparent  stresses  and  their  derivatives  at  various 
mesh  points.  The  procedure  for  calculating  these 
derivatives  is  the  same  as  that  described  for 
subroutine  ROEN  2.  This  subroutine  is  called  by 
specifying  input  variable  NWAT  =  1.0  . 

The  logic  of  the  above  subroutines  and  their  interaction 
with  the  main  program  is  given  in  Figure  25. 

6 . 3  Computation  of  Mean  Flow  Quantities  and  the  Apparent 

Stresses  from  the  Pxpsri mental  Data 

The  program  named  EXPTL  FORT  follows  the  following 
sequence  of  steps  for  estimating  mean  flow  quantities  and 
the  apparent  stresses. 

(i)  Read  experimental  data  on  the  instantaneous 
values  of  Mach  numbers  Mr2,  M©2,  Mz2#  at  static 
pressure  at  different  radial  locations.  Read 
other  input  data  such  as  inlet  total  temperature 
and  pressure. 

(ii)  First  of  all  the  total  temperature  ratio  is 
determined  by  using  the  Euler  Turbine  equation 
in  the  following  form: 
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where  mt  is  the  rotor  tip  Mach  number  and 
Mz 

(iii)  With  Tt^  known,  Tt2  and  TSj  can  be  estimated. 
Instantaneous  values  of  the  velocity  of  sound 
and  the  flow  velocities  are  now  estimated  at 
different  radial  locations. 

(iv)  Mean  velocities  and  pressure  are  determined  by 
using  the  relation 

t-  V  Mi 
V  =  * 

(v)  Once  the  mean  velocities  have  been  estimated, 

the  fluctuating  components  of  velocities  can  be 
estimated  1^"  =  \}  -  V 


SECTION  VII 


RESULTS  OP  NUMERICAL  CALCULATIONS  AND  DISCUSSION 

The  mean  flow  field  of  the  MIT  Blowdown  Rotor  was 
computed  by  using  the  modified  Streamline  Curvature  Computer 
Program,  described  in  the  preceding  chapter.  The  purpose  of 
this  exercise  was  to  verify  the  validity  of  flow  models  for 
Apparent  stresses.  Mean  Rothalpy  variation  along  the  stream¬ 
line,  and  Apparent  Entropy  Production,  plow  models  for  these 
phenomena  which  are  associated  with  pitchwise  flow  fluctua¬ 
tions,  were  formulated  in  Section  III. 

First  of  all,  the  mean  flow  field  of  the  MIT  Blowdown 
rotor  was  computed  by  only  using  the  cascade  losses.  This 
was  to  confirm  that  the  original  streamline  Curvature  program 
is  able  to  accurately  predict  the  design  flow  field  (free 
vortex) .  Computations  were  next  carried  out  to  study  the 
effect  of  the  above  three  phenomena  on  the  mean  flow  field. 
Here  the  effect  of  individual  phenomenon  and  the  combined 
effect  of  these  phenomena  are  studied,  computations  were 
also  carried  out  by  superimposing  the  effect  of  shock  losses 
on  the  above  calculation.  Finally,  the  flow  field  was  cal¬ 
culated  by  using  experimentally  estimated  losses. 

The  apparent  stress  distribution  given  in  Figures  9 

through  20  were  used  as  input  for  estimating  II  ,7*2^  1 

9m  *pp 

and  the  Apparent  Stress  terms  appearing  in  the  equation  of 
motion  (50) . 


The  negative  entropy  rise  through  the  rotor  as  seen 
around  the  midspan  region  (Figure  22) ,  made  the  problem  of 
specifying  the  experimental  entropy  distribution  a  rather 
difficult  task.  In  the  absence  of  any  reasonable  explanation 
this  negative  entropy  rise  was  attributed  to  a  zero  shift  in 
the  pressure  transducer.  Therefore,  for  the  purpose  of 
specifying  the  experimental  losses,  the  curve  of  Figure  22 
was  shifted  up  so  as  to  bring  the  most  negative  value  to  zero. 
This  is  shown  in  Figure  26. 

7.1  Computational  Mesh  and  Basic  Input  Data 

The  computational  mesh  for  application  to  the  MIT  Blow¬ 
down  compressor  was  formed  by  seven  nearly  radial  fixed 
stations  intersecting  eleven  meridional  streamlines.  Two 
stations  were  located  ahead  of  the  rotor,  one  each  at  the 
rotor  leading  and  trailing  edges,  two  within  the  blade  row, 
and  one  at  a  distance  of  1.0  chord  downstream  of  the  trailing 
edge.  The  hub  and  the  casing  form  two  fixed  streamlines. 

The  positions  of  the  remaining  nine  streamlines  are  fixed 
during  any  computational  pass  but  are  readjusted  at  the  end 
of  each  pass.  Figure  27  gives  the  details  of  the  computation¬ 
al  mesh.  It  may  be  recalled  that  the  time  resolved  experi¬ 
mental  data  was  obtained  at  0.1  chord  (Port  5)  and  1.0  chord 
(Port  6)  downstream  of  the  rotor.  For  this  computational 
study  the  port  5  location  was  taken  to  represent  the  blade 
trailing  edge  projection  on  the  meridional  plane. 
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The  basic  input  data  for  the  anlysis  problem  under 
investigation  consisted  of  blade  exit  relative  flow  angles, 
geometrical  features,  mass  flow  rate,  rpm,  and  the  compressor 
inlet  conditions.  For  all  the  computational  studies,  the 
experimentally  determined  relative  exit  flow  angles  ( y3  ) 
were  used.  Alternately,  the  flow  angles  can  be  specified  by 
using  the  deviation  angles  based  on  the  cascade  data.  Al¬ 
though  this  does  not  alter  the  basic  character  of  the  flow 
field,  the  experimental  flow  angles  were  used  to  get 
better  accuracy.  Relative  flow  angles  for  stations  within 
the  blade  row  were  calculated  by  assuming  a  sinusoidal  dis¬ 
tribution  of  cLC'l'Vq)  along  the  axial  chord,  it  may  be 
mentioned  here  that  any  alternate  load  distribution  within 
the  blade  row  will  have  negligible  effect  on  the  exit  flow 
field,  provided  the  load  distribution  at  the  rotor  exit 
remains  unchanged.  This  is  because  a  change  within  the  blade 
row  can  be  communicated  to  the  blade  exit  location  only 
through  the  change  of  the  slope  and  the  curvature  of  stream¬ 
lines  which  in  turn  are  primarily  controlled  by  the  load 
distribution  at  the  rotor  exit.  Other  basic  input  data  such 
as  the  mass  flow  rate  and  rpm  were  reduced  to  the  standard 
atmospheric  conditions  at  the  compressor  inlet  with  air  as 
the  fluid.  Details  of  the  basic  input  data  and  apparent 
stress  data  are  given  in  Figures  28a  and  28b. 
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7.2  Mean  Flow  Field 


7.2.1  Computed  Flow  with  Cascade  Losses 

Mean  pitchwise,  axial  and  radial  velocity  distributions 
at  the  rotor  exit  (Station  6)  and  at  station  7,  as  computed 
by  using  cascade  losses,  are  shown  in  Figures  29  through  32. 

It  is  seen  that  the  computed  flow  field  is  significantly 
different  from  that  obtained  experimentally.  The  pitchwise 
velocity  has  an  almost  constant  negative  slope  from 

the  root  to  the  tip  and  has  completely  failed  to  predict  the 
double  bucket  shaped  distribution  found  experimentally,  in 
fact,  the  steep  downward  trend  of  near  the  casing  would 

cause  flow  instability  and  is  therefore  not  expected  to  exist 
in  the  actual  turbomachine. 

The  computed  distribution  of  mean  axial  velocity  does  not 
show  a  steep  fall  near  the  casing,  which  is  an  important 
feature  of  the  experimental  14  distribution. 

The  radial  velocity  distribution  which  is  mainly  con¬ 
trolled  by  the  large  slope  of  the  hub,  is  found  to  decrease 
steeply  from  a  large  value  at  the  hub  to  zero  at  the  casing 
in  conformity  with  the  casing  boundary  condition.  A  reversal 
of  slope  at  two  locations  as  observed  experimentally,  has  not 
been  predicted  by  this  computation. 

It  may,  however,  be  noted  that  the  results  of  this  com¬ 
putation  show  a  close  resemblance  to  the  flow  field  predicted 
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by  the  three-dimensional  inviscid  computational  scheme  of 
Reference  (2) .  The  predicted  flow  field  also  gives  a  good 
agreement  with  the  design  flow  field  of  free  vortex  type 
thus  showing  that  the  basic  program  of  Reference  (1)  is 
reliable  for  the  prediction  of  mean  flow  field.  Comparison 
with  the  3-D  inviscid  computation  also  shows  that  3-D  invis¬ 
cid  effects  per  se  do  not  alter  the  mean  flow  field  of  this 
rotor  in  a  major  way.  The  principal  3-D  effects  are  due  to 
viscous  interactions,  vfaich  are  of  course  strongly  controlled 
by  the  inviscid  3-D  flow. 


7.2.2  Computed  Flow  with  Cascade  Losses  and 
Apparent  stresses 


Figures  33  through  36  give  the  mean  flow  field  com¬ 
puted  with  cascade  losses  and  Apparent  Stresses.  Variations 
of  mean  Rothalpy  and  the  production  of  Apparent  Entropy  along 
the  meridional  streamline  are  not  included  in  this  computa¬ 
tion.  It  is  seen  from  Figure  33  that  the  mean  pitchwise 
velocity  distribution  at  the  rotor  exit  is  similar  to  that 
obtained  without  Apparent  stresses,  a  slight  depression  at 
r/^  =  0.87  and  a  slight  increase  in  the  negative  slope  near 
the  casing  indicate  a  trend  reverse  of  that  observed  experi¬ 
mentally.  Although  the  direct  effect  of  Apparent  stresses 

is  negligible  on  the  mean  flow  field,  it  is  noted  from 

__2 

Figure  37  (app.  stress/  Vfe  vs  r/  that  the  order  of 
magnitude  of  Apparent  Stresses  is  comparable  to  that  of  the 


.  A  comparison  of  Figure  37  with 


centrifugal  force 

- —  T. 

Figure  12a  ( vf  't/j^ 

all  the  contribution  to  the  Apparent  stress  terms  of  equa¬ 


)  shows  that  almost 


tion  (50)  comes  from  the  distribution  of  the  normal  radial 
Apparent  Stress 


7.2.3  Computed  Flow  with  Cascade  Losses  and 
Apparent  Entropy  Production 


Before  computing  the  mean  velocity  with  cascade  losses 
and  Apparent  Entropy,  the  distributions  of  7"^  )  at 
stations  4  and  5  which  are  located  within  the  blade  row, 
were  estimated.  Apparent  Entropy  rise  across  the  rotor  is 
then  computed.  These  distributions  are  given  in  Figures  38 
through  40.  It  is  seen  from  these  figures  that  the  appa¬ 
rent  entropy  production  is  mainly  controlled  by  the  term 
j&%).  Since  the  absolute  value  of  apparent  stress  and 
the  mean  velocity  gradients  are  small,  the  apparent  dissipa¬ 
tion  term  S  in  equation  (32)  is  small. 


The  computed  mean  flow  field  with  cascade  losses  and 
Apparent  Entropy  is  given  in  Figures  41  through  44. 

It  is  seen  from  Figure  41  that  in  the  midspan  region, 
at  the  rotor  exit  shows  a  trend  similar  to  the  experimental 
results.  It  is,  however,  noted  that  the  extent  of  change 
introduced  by  Apparent  Entropy  in  the  computed  flow  is 
considerably  less  than  observed  in  the  experimental  distribution. 


7.2.4  Computed  Flow  with  Cascade  Losses  and 


Rothalpy  Variation  Along  the  Streamline 


*3/ 

Distributions  of  -  at  stations  4  and  5,  and  mean 

Rothalpy  variation  across  the  rotor  are  first  computed  and 
these  are  given  in  Figures  38,  39,  and  45.  It  is  seen 
from  these  figures  that  mean  Rothalpy  decreases  steeply  to¬ 
wards  the  hub  and  the  casing,  and  it  increases  slightly  in  the 
midspan  region.  This  distribution  is  due  to  the  spatial 


transfer  of  mean  flow  kinetic  energy  from  the  hub  and  the 
casing  towards  the  midspan  region.  Besides  towards  the 
casing,  conversion  of  mean  flow  kinetic  energy  to  the  energy 
of  fluctuations  becomes  quite  significant. 


The  computed  mean  flow  field  with  cascade  losses  and 
Rothalpy  variation  is  given  in  Figures  46  through  50. 

The  shape  of  the  curve  at  the  rotor  exit  (station  6) 

shows  good  agreement  with  the  experimental  distribution  up  to 
r/rt=  0.9.  It  is  further  seen  that  a  good  agreement  exists 
between  the  predicted  and  the  experimental  distributions  of 
and  at  station  7.  Thus  the  effect  of  mean 

Rothalpy  variation  along  the  streamline  seems  to  play  an 
important  role  in  changing  the  mean  velocity  field. 

7.2.5  Computed  Flow  with  Cascade  Losses,  Rothalpy 
Variation,  and  Apparent  Entropy  production 
Along  the  streamline 

rn  order  to  visualize  the  combined  effect  of  Apparent 
Entropy  production  and  mean  Rothalpy  variation,  equivalent 
entropy  change  across  the  rotor  defined  as 
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^  ~  ) 

is  computed.  This  is  plotted  in  Figure  51.  It  is  interesting 
to  note  that  a  close  resemblance  exists  between  the  equivalent 
entropy  rise  and  the  experimental  entropy  rise. 

Figures  52  to  56  give  the  computed  mean  velocity 
field  at  station  6  and  7.  The  combined  effects  of  Rothalpy 
Variation  and  Apparent  Entropy  production  is  found  to  change 
the  free  vortex  velocity  distribution  to  a  distribution  which 
closely  resembles  the  experimental  distribution.  The  obser¬ 
ved  difference  in  the  absolute  values  of  Va  at  station  6, 
is  attributed  to  the  lower  experimental  values  resulting  from 
a  zero  shift  in  the  pressure  transducer.  Therefore,  the  shape 
of  the  curves  should  only  be  compared  to  verify  the  validity 
of  the  flow  model. 

7.2.6  Commuted  Flow  with  cascade  Losses,  Rothalpy 
Variation,  Apparent  Entropy  Production,  and 
Apparent  Stresses 

Finally,  the  flow  field  was  computed  with  all  of  the 
above  mentioned  phenomena,  and  the  results  of  numerical 
calculations  are  given  in  Figures  57  through  61.  It  is 
observed  that  the  effect  of  adding  Apparent  Stress  to  Roth¬ 
alpy  variation  and  Apparent  Entropy  production  is  small,  and 
is  localized  near  the  casing.  In  fact  the  positive  slope  of 
v9  which  was  obtained  by  using  Rothalpy  variation  and 
Entropy  production,  no  longer  exists  here,  several  factors 


could  be  responsible  for  the  disagreement  between  the  experi¬ 
mental  and  the  predicted  results  near  the  casing.  Error  in 
the  experimental  data  near  the  casing,  assumption  of  linear 
distribution  of  apparent  stresses  along  the  meridional  stream¬ 
line,  and  non- inclusion  of  shock  losses  are  some  of  the  im¬ 
portant  factors.  The  effect  of  superimposing  shock  losses 
on  the  above  computation  is  given  in  Figure  62.  This  is 
seen  to  make  the  slope  of  'V#  positive  near  the  casing. 

7.2.7  Computed  Flow  with  Experimental  Losses 

Mean  flow  field  computed  by  using  the  entropy  distribu¬ 
tion  of  Figure  26  is  given  in  Figures  63  through  68. 

It  is  observed  that  Vq  predicted  at  the  rotor  exit  has  an 
excellent  agreement  with  the  experimental  results. 
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SECTION  VIII 


CONCLUSIONS 

The  primary  aim  of  this  research  has  been  to  improve 
the  existing  turbomachine  flow  model  based  on  the  combina¬ 
tion  of  ax i symmetric  through  flow  and  2-d  cascade  flow,  so 
that  it  can  be  applied  to  highly  loaded  systems  for  which 
the  flow  field  computed  by  the  existing  flow  model,  has 
shown  poor  agreement  with  the  experimental  results.  In  order 
to  achieve  this  goal,  the  existing  flow  model  has  been  modi¬ 
fied  to  account  for  blade  to  blade  flow  variations.  The 
major  support  in  verifying  the  improved  flow  model  came  from 
the  time  resolved  experimental  data  of  the  MIT  Blowdown 
Compressor  facility.  A  critical  analysis  of  this  data  con¬ 
siderably  helped  in  understanding  the  3-D  viscous  effects 
inherent  in  a  highly  loaded  axial  flow  compressor.  The  com¬ 
puted  flow  field  obtained  by  modifying  the  streamline  Curva¬ 
ture  program  of  Reference  (1)  proved  useful  in  consolidating 
many  of  the  theoretical  concepts  and  led  to  several  important 
conclusions.  A  general  disucssion  on  the  highlights  of  this 
work  and  the  conclusions  arising  out  of  it  are  summarized  here. 

8.1  In  Sections  I  and  ill  of  this  report  it  is  shown  that 
the  effect  of  blade  to  blade  flow  variations  on  the  mean 
(circumferentially  averaged)  flow  can  be  represented  by  the 
following  phenomena: 
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(i)  Apparent  stresses  resulting  from  the  pitchwise 
averaging  of  the  momentum  equation. 

(ii)  Mean  Rothalpy  variation  along  the  streamline  due 
to  energy  transfer  across  the  stream  surfaces  by 
by  the  apparent  stresses,  and  by  conversion  of 
mean  flow  kinetic  energy  to  the  energy  of  fluc¬ 
tuations  . 

(iii)  Apparent  Entropy  Production  due  to  dissipation 
associated  with  apparent  stresses,  and  due  to  the 
production  of  fluctuation  energy. 

These  phenomena  significantly  influence  the  spatial 
distribution  of  mean  momentum  and  energy  in  an  axial  flow 
compressor. 

8.2  in  Section  V  a  critical  study  of  the  time  resolved  ex¬ 
perimental  data  of  the  MIT  Blowdown  Compressor  is  carried  out. 
Some  of  the  important  findings  of  this  study  are: 

(i)  A  strong  disagreement  exists  between  the  mean 
experimental  flow  field  and  the  design  flow  field. 

(ii)  Apparent  stresses  resulting  from  the  pitchwise 
flow  variations  are  of  significant  magnitude. 

T* 

in  particular,  the  radial  derivative  of 

at  the  rotor  exit  is  found  to  be  of  comparable 

magnitude  to  the  mean  specific  centrifugal  force 
_ 2 

Vo  .  Since  the  normal  radial  apparent  stress 
represents  the  spanwise  flows  in  the  blade  boun¬ 
dary  layer,  the  above  observation  indicated  that 
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for  a  highly  loaded  system,  the  radial  transport 
of  boundary  layer  and  wakes  is  an  important  effect. 

(iii)  The  radial  gradient  of  the  mean  entropy  rise  near 
the  hub  and  the  casing  is  estimated  to  be  several 
times  larger  than  that  predicted  by  two-dimension¬ 
al  cascade  data. 

8.3  a  detailed  study  of  the  effect  of  phenomena  listed  in 
Section  8.1  has  been  carried  out  by  using  the  flow  models 
developed  in  Section  III.  A  discussion  on  the  flow  field 
computed  with  the  aid  of  modified  streamline  curvature  compu¬ 
ter  program  (Section  VI)  is  presented  in  Section  VII,  and 
the  conclusions  arising  from  this  discussion  are  summarized 
as  follows: 

(i)  An  axisymmetric  computation  with  entropy  distri¬ 
bution  based  on  the  cascade  data  shows  an  excel¬ 
lent  agreement  between  the  predicted  flow  and  the 
design  flow  and  also  agrees  well  with  the  results 

of  a  3-D  inviscid  computation.  It  is  therefore  con¬ 
cluded  that  3-D  inviscid  effects  such  as  the  blade 
passage  shock  and  the  distortion  of  stream  surfaces 
due  to  the  inviscid  pressure  field  of  the  blades, 
do  not  produce  significant  changes  in  the  mean  flow 
for  this  rotor.  Hence,  the  deviation  from  the  de¬ 
sign  flow  of  the  experimental  results  is  attributed 
to  3-D  viscous  effects  which,  however,  are  indirectly 
controlled  by  the  inviscid  blade-to-blade  flow. 

(ii)  The  introduction  of  apparent  stresses  into  the 
mean  flow  1-direction  (nearly  radial)  momentum 
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equation  does  not  have  a  large  effect  on  the 

mean  flow  field,  it  is,  however,  found  that  the 

major  contribution  to  the  apparent  stress  terms 

in  the  above  moment  vim  equation  comes  from  the 

“7* 

radial  derivative  of  .  in  fact  a  knowledge 

of  this  derivative  which  represents  the  distribu¬ 
tion  of  spanwise  flow  in  the  blade  boundary  layers 
and  wakes,  will  probably  be  adequate  to  replace 
all  the  apparent  stress  terms. 

(iii)  Mean  Rothalpy  variation  and  the  production  of 
Apparent  Entropy  along  the  streamline  play  the 
major  roles  in  reshaping  the  velocity  field  of  a 
highly  loaded  axial  flow  compressor.  The  flow 
field  computed  by  accounting  for  these  effects 
shows  a  distinct  resemblance  to  the  experimental 
flow.  The  largest  contribution  in  reshaping  the 
velocity  field  comes  from  the  variation  of  mean 
Rothalpy  along  the  streamline.  The  effect  of 
Apparent  Entropy  production,  however,  is  consider¬ 
ably  larger  than  that  due  to  Apparent  Stresses  and 
hence,  it  cannot  be  ignored,  it  is  thus  concluded, 
if  the  mean  Rothalpy  variation  and  the  Apparent 
Entropy  production  along  the  streamline  are  taken 

into  account,  the  axisymmetric  computational 
technique  will  be  able  to  predict  the  flow  field 
of  a  highly  loaded  axial  flow  compressor  with 
reasonable  accuracy. 
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APPENDIX  A 


DERIVATION  OF  RADIAL,  TANGENTIAL  AND  AXIAL 
EQUATIONS  OF  MOTION 

If  we  choose  a  frame  of  reference  fixed  to  the  rotor 
blades,  unsteadiness  in  the  flow  exists  due  to  turbulent 
fluctuations  only.  Equations  of  mass  and  momentum  conser¬ 
vations  for  unsteady  flow  are: 

Mass  Conservation 


+  .L  Wt)  /  ^(OVe)  'ZCfWz) 


T>t  '  Z 

momentum  Conservation 
r-direction 


=  0 


9-direction 


--L2L  +  ?** 

f  &L 


lr  +  wjW:  +aJi  Yh 


z-direction 


=  -i-3^  -f-  F&v- 

?  7& 


•jr  +  m  ^ _ F 

at  ^  ^ —  —  ~  — *r 


f  2* 


(A-2 ) 


(A-3) 


( A-4 ) 


in  order  to  make  these  equations  time  independent,  a  time 
average  of  these  equations  will  be  carried  out,  which  is  de¬ 
fined  by  the  following  integral:  j. 

A  >  —  I A*± 


(A- 5) 
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✓A  =■  +  A 


(A-5a) 


Time  averaged  mass  conservation  equation  becomes 

_/  jH<^>  —  2.^We>  _ 

't  9<fc  ^90  ^  7£.  ° 


(A-6) 


Here  is  the  time  averaged  mass  flow  defined 

J 

=  c^/m  —  f(tw)cte 
T-*°«  T  i 

•  V  f  •*.  n  ' 


fVV  +  ffw)' 


(A- 7) 


(A- 7  a) 


For  derivation  of  the  time  averaged  momentum  equations, 
the  mass  conservation  equation  is  first  combined  with  the 
momentum  conservation  equation  followed  by  time  averaging  of 
the  resulting  equation.  The  time  averaged  radial  momentum 
equation  becomes 

2£y*->  +  _<Wo>sl 

dt  ^  DO 

-  ^  Jl2  _  ^  we>  _  _  J__  f+v? 

*  ^  <?;> 


+  Mt  >*  _  Wp)Ve/^> 


(A- 8) 


104 


In  order  to  account  for  the  blade  to  balde  flow  varia¬ 
tions,  pitchwise  averaging  of  this  equation  is  carried  out 
as  per  the  procedure  outlined  in  Appendix  B.  Denoting  pitch- 
wise  average  of  the  time  averaged  quantities  by  <  ^  ,  and 

deviation  from  pitchwise  average  by  ?  ,  the  mean  mass 

conservation  and  the  radial  momentum  equations  are  derived 
as  follows: 


Mass  Conservation  Equation  (Pitchwise  Averaged) 

Following  Appendix* B,  pitchwise  average  of  equation  (A-ll) 
is  represented  as: 


,  Qj  Or  Or  Oj 

+  J -«r d&  +  Jm  l*‘=o 

I  II  III  IV 

Individual  terms  can  now  be  evaluated. 

From  equation  (BIO)  of  Appendix  B,  Term  I  becomes: 


Term  I  - 


(f<V^  -  (< t£») ? 

)  /  2> 

\  2.  /  > 


(A-13) 


(A-14 ) 


Assuming  no  slip  condition  on  the  blade  surfaces, 
(A14)  becomes 


Term  i 


s  ^  vUr 

'M. 


1 

*  dt 


equation 

( A- 1 5 ) 
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Here  W  (or  (  f W  )'  )  represents  temporal  fluctuations 
about  the  time  averaged  value  represented  by  iv>  (or  <£5  yy;>)  • 
The  third  term  in  the  right  hand  side  of  equation  (A-8)  gives 
the  contribution  of  turbulent  stresses. 


An  alternate  method  of  deriving  time  averaged  equations 
of  motion  is  to  separate  density  fluctuations  from  the  mass 
flow  fluctuations  used  in  the  above  derivation.  The  time 
averaged  mass  conservation  and  radial  momentum  equations 
obtaining  by  using  this  procedure  are: 

Mass  conservation 


»  V  - -  ►  / 


J  o  V 

3t  't  o<o 

+ U  2.  -3  j 

l*71  J 


=  o 


(A-9) 


Radial  Equation  of  Motion 

^ «,>  2<"“- V  ^  j2 

^  ^30  3^ 

<2*#  =  — /,  « w 

+  +  <ci'w£?’U 

*  ^  {—+<77^Jy 

-l-p,  >«*> 

+2<Mzn 1  ] 


(A-10) 
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The  last  term  in  the  right  hand  side  of  equation  (A-10) 
is  the  apparent  stress  contribution  from  the  time  averaged 
mass  conservation  equation. 


Both  the  equations  (A-8)  and  (A'-10)  require  information 
on  the  instantaneous  value  of  density.  But  since  the  exist¬ 
ing  experimental  methods  cannot  provide  reliable  data  on 
density,  it  .was  considered  desirable  to  derive  the  mean  flow 
equations  for  the  case  of  small  density  fluctuations.  The 
time  averaged  equations  of  mass  conservation  and  radial  mo¬ 
mentum  now  become: 

Mass  Conservation  Equation 


\  '  °  (A- 11) 


Radial  Momentum  Equation 

t  90  2  3  * 

-2<«*>SL  -'tSt-  ^  7 

_  /^Kh4  w) _/  [wLx) 

I'M,  + 


*  90 


9* 
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Similarly  treatment  gives 


Term  II 


■4WtM 


Term  III  =  Q 


Term  IV 


_  wz>; 

H>z 


7  U 


>  (A-16) 


Thus  the  pitchwise  averaged  mass  conservation  equation 
becomes : 


4 rnv>  c»i[>  "PA  ,  -£v/b>  pA 

+~  +—  S  = 

Multiplying  equation  (a-17)  by  ,  gives 

*  =  o 


Equation  (A-18)  represents  the  axisymmetric  form  of  the 
mass  conservation  equation  in  which  the  blockage  due  to  the 
blade  thickness  has  been  uniformly  distributed. 


in  deriving  the  above  equation  and  also  for  the  subse¬ 
quent  derivations,  it  is  assumed  that  the  pitchwise  variation 
of  density  is  small. 

Radial  Momentum  Equation  (Pitchwise  Averaged) 

Multiplying  equation  (A-ll)  by^'Wt>  and  adding  it  to 
equation  (A-12),  gives 
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■3.  r  ^  ^  *l 

-2<wd>ji?_fsf't- 1  2*p>  /-  r^/  i»  , 


+4£?*‘*’+£"te’i  <<*-<<> 


(A-19 ) 


Pitchwise  averaging  of  individual  terms  can  now  be  carried 
out  using  equations  (B-ll) ,  (B— 16 ) ,  and  (B— 1 7 ) 


'2kvj%>ty _  '3^>f  -JX  .  "3 ^OT^>f  -jv 

”*•  -  *  +-j~w^  +  —  If 


(A-19a) 
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(  ^W/t> 

ru 
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^  't  <fc 

(A-22a) 

-  r 

( A-  2  3 ) 

109 

(  _ 

where  Fri  is  the  inviscid  blade  force  component  in  the  radial 
direction  and  is  equal  to 

2  71  Vn 

We  can  thus  write  the  pitch  averaged  form  of  equation  (A-19) 

^  +2j£y/*>’,)+  2b  ^  («*>  ? ) 


_ — -  y  > 

>  >t 


^  *32  ^ 


>  5i 


+?<K  - 


IK  -  ^We>f 

*  *  ~ 


=  -  ^^+^1  +f<*W>  +fft. 


(A- 2 4 ) 


Here  Frt  is  the  contribution  of  pitchwise  averaged  turbulent 
stresses,  given  by 


rn  =  -il -<60 


V. 


(A- 2  5 ] 


Multiplying  equation  (A-24)  by  /•  and  after  simple  algeb] 


man ipu 1 at  ion s ,  we  get 


^  J^<r yvx>^^>+  <Wq>_2SL  <wq>  -S)X 


+  <W-t> 


2_(><vvz>^ 


(A- 26) 
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+  [  <"£>>)  +  >^we>f 

✓ 

=  a/-  25>+fA<  +r<^;>  +rSi  j 

Usirg  mass  conservation  equation (A-18)  and  dividing  both 
sided  by  ,  we  get  final  form  of  pitch  averaged  radial 

momentum  equation 


32 


^  -  2S1<Wq> 


~  =  -  4-  — ^ •  ±-  /V^<vVt>$  3  (Vwt> <wa>fj_ A  <w^l 

T  ^  f>(3t  32 


+  <  fikv>  -H  Rt 


( A- 2  7 ) 


This  is  the  final  form  of  radial  momentum  equation 
which  has  first  been  time  averaged  to  make  it  time  indepen¬ 
dent,  and  then  pitchwise  averaged  to  make  it  9- independent . 

Equation  (A-27)  contains  two  types  of  velocity  correla¬ 
tions,  one  accounting  for  turbulent  fluctuations  (temporal) , 
and  the  other  accounting  for  the  spatial  fluctuations  of  the 
time  averated  quantities.  For  a  highly  loaded  system  where 
blade  to  blade  flow  variations  are  large,  the  contribution 
from  turbulent  fluctuations  is  small  as  compared  to  that  due 
to  the  pitchwise  fluctuations. 


Ill 


The  above  procedure  can  also  be  applied  to  the  flow 
fluctuations  outside  the  blade  row  in  the  absolute  frame 
of  reference  to  account  for  pitchwise  flow  variations. 

This  is  done  by  first  taking  an  ensemble  average  of  the  un¬ 
steady  equation  of  motion  and  then  taking  time  average  of 
ensemble  averaged  quantities. 


To  simplify  equation  (A-27),  we  can  either  neglect  tur¬ 
bulent  stresses  as  compared  to  the  large  scale  blade  to  blade 
fluctuations,  or  a  slightly  different  averaging  procedure 
can  be  followed,  it  is  possible  to  directly  carry  out  pitch- 
wise  averaging  of  unsteady  equations  of  motion  by  assuming 
to  be  small.  This  procedure  will  give  the  following 


'it 

radial-direction  equation  of  motion  in  the  apparent 

£°mat'  Yk  ^  +  -  ZWg  J1  - 

Ot  'Sz  't  _ _ 

=  _ 1 2A  ±.f± 9 [h ** »£?)_ I 

?  3t  'L  ) 


stress 


(A-28) 


+ 

Corresponding  equations  in  9  and  z  directions  are 


9-direction 


wt  ^1°  +  wz  +  w&Wt  +  2-&  Wi 

_ _ -  -fc.  _  _ 

-  _  l  f WhIms*)  Ki'lW  +  Fa,  +  Fov  (A-29) 
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z-direct ion  _  _ 

Wt 
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APPENDIX  B 


PITCHWISE  AVERAGING 


Pitchwise  average  of  any  flow  property  a  within  the 
blade  row  is  defined  as  a 


A  = 


A  ctQ- 


(B-l) 


Here  the  integration  has  been  carried  out  from  the  pressure 
side  of  one  blade  to  the  suction  side  of  the  adjacent  blade. 
This  is  consistent  with  the  direction  of  rotation  for  an 
axial  compressor  rotor.  For  a  turbine  rotor  the  integration 
would  be  carried  out  from  the  suction  side  of  one  blade  to 
the  pressure  side  of  the  adjacent  blade. 

Following  Smith**,  the  mean  blade  surface  and  the  blade 

lean  angle  U  is  defined  in  Figure  (3.1).  if  the  blade 

blockage  factor  is  now  defined  as  A-  *  equation 

•s 

(Bl)  becomes 


A  = 


AciQ- 


(B-2) 


Here  (  /—  A  )  represents  the  blockage  to  the  flow  due  to  the 
blade  thickness. 

in  order  to  carry  out  pitchwise  average  of  derivative 
the  flow  property  we  adopt  Liebnitz  Rule  for  expressing  the 
derivative  of  an  integral. 

Thus 
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Pitchwise  Averaging  of 


>/3 


=5 


Equation  (B-3)  will  now  be  used  to  carry  out  pitchwise 


averaging  of  2k 

h'c 


2k  = 
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Term  I 
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For  evaluating  Term  n,  ©s  and  0p  are  expressed  as 
follows 


©,  =  Sw  -  Z£  0->~J  +  2k 

'Y  /V 


(B-6) 
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ef=  =  -  js  o-y) 


here  ©m  is  the  angular  location  of  the  mean  blade  surface. 
From  Figure  (3.1)  the  blade  lean  angle  is  given  by 


tun  s  -  t 


(B-8) 


Using  equations  (B-6) ,  (B-7)  and  (B-8) ,  Term  II  becomes 


Term  II 


=  -  O--^ 

=  -  Qi±b)  _L  It  .  (t>* -f>P) 


2  >  2t  ^ 

Therefore  equation  (B-4)  becomes 
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(B-10) 


The  term  kJjL-L£-  defines  the  radial 

2-TV>/N 

component  of  the  inviscid  blade  force. 

The  term  yb  —  ~T~  will  be  small  unless 

is  very  large.  For  the  practical  blade  geometries 
is  small;  fr-  ^ JA  is  also  small.  Hence,  the  last 
term  of  equation  (B-10)  can  be  neglected,  giving 


will  be  small  unless 

^ _ _  3* 


*2t  _  'djb 


^  = 
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Here,  terms  linear  in  the  fluctuating  quantities  have 
not  been  retained  since  these  will  vanish  on  pitchwise 
averaging. 

Following  the  pitchwise  averaging  procedure  for 
and  noting  that  velocities  at  the  blade  surface  will  vanish 
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(no  slip  condition),  we  get 


APPENDIX  C 


MEAN  FLOW  ENERGY  EOUATION 
(Frame  of  Reference  Fixed  to  the  Blades) 


Assumptions: 

(i)  Adiabatic  flow 

(ii)  Blade  to  blade  variations  in  density  and 
internal  energy  are  small 

(iii)  Work  done  by  viscous  stresses « work  done  by 
apparent  stresses 

(iv)  Mean  flow  is  invarient  in  time 
1.  Rate  of  Gain  of  internal  Energy  per  Unit  volume 

=  2(<eyk)  f£Wt  j_  2(<?we)  2_(!£wi) 

t  ^95  9z 


(c-i) 


Where  E  is  the  internal  energy  per  unit  mass, 
pitchwise  average  of  equation  (C-l)  gives: 

Rate  of  gain  of  mean  internal  Energy/vol 

_  5>jr 

where  J>_  =  9.  y-  Q 

3>t 

2 .  Rate  of  Gain  of  Kinetic  Energy  per  Unit  Volume 

_2< +±2 (r%4)^(fw2wX) 


Taking 


(C-2) 


(C-3 ) 
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Where 


( C— 5 ) 


w*=  w£ 


Using  continuity  equation  V»(?iv}=0  gives: 
Rate  of  Gain  of  K.E/vol 

_  ?  vwj  2K/y+  ?  iva  J  5  w  9^; 
^  T  'Se  -32 

Taking  pitchwise  average  gives: 


Rate  of  Gain  of  mean  K.E/vol 


(C-7) 


(C-8) 


3.  Rate  of  work  Input  from  Pressure  Forces  per  Unit  Volume 


The  second  term  in  the  right  hand  side  of  equation  (C-10) 
can  be  written  as  V»  W  and  using  continuity  equation, 
we  have 
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(C— 12) 


4.  We  can  find  the  work  done  per  unit  volume  by  the  normal 
and  shear  stresses.  Here  we  will  assume  that  the  apparent 
stresses  introduced  due  to  the  pitchwise  averaging  of  equa¬ 
tions  of  motion  dominate  and  the  viscous  stresses  are  small. 

Work  done  by  apparent  normal  stresses  per  unit  volume 


= -j/*^**£ §  j 

Here  the  terms  containing  the  derivative  of  density  are  neglected. 


■ML 


(C-15) 


Similarly,  work  done  by  apparent  shear  stresses  per  unit 


If  and  —  represent  the  apparent  stress  terms 

appearing  in  the  radial,  pitchwise  and  axial  momentum  equa¬ 
tions  respectively,  and  if  j?  denotes  the  production  or 
apparent  dissipation  term,  the  rate  of  work  done  by  apparent 
stresses  per  unit  volume  can  be  written  as  : 


Rate  (work  done  by  apparent  stresses) 

=  j  +  (C-17) 

where  $  app  =  -  {  2^ 

v*‘t 


5.  Rate  of  Work  Done  by  centrifugal  Force 

As  we  are  using  the  relative  frame  of  reference,  it  is 
necessary  to  account  for  the  work  done  by  the  centrifugal 
force.  it  may  be  mentioned  here  that  no  work  is  done  by  the 
Coriolis  forces  since  both  the  radial  and  the  tangential 
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components  of  coriolis  forces  are  associated  with  velocities 
normal  to  them. 

Centrifugal  force  only  has  a  radial  component  and  there¬ 
fore  work  done  by  it  per  unit  volume  z 

(C-18) 

Pitchwise  averaging  gives  the  mean  rate  of  work  done  by 
the  centrifugal  force 

=  Wfc  'tsi1-  = 


3>t 


(C-19 ) 


The  energy  equation  can  now  be  written  by  equating  the 
rate  of  gain  of  mean  internal  and  kinetic  energies  to  the 
mean  work  input  rate  from  the  pressure  and  centrifugal  forces, 
and  the  apparent  stresses.  This  gives,  after  rearranging 

%  _  pILC^A) 

x>t  X>t  'TU.  a  STT  j 


3>t  ? 


*3>t 


=  +  +  waJzj-  $ 

e  f  3>'w,l/2  d'C wwD  l  _  ^Cw'%] 

U  a>t  a>t  /  3* 


(C- 20) 


-f* 


Now,  consistent  with  the  discussion  presented  in  Section 


3.3,  mean  rothalpy  is  defined  as 

T  =  E  -h  JL  + 

Z  2-2. 


(C-21) 
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Equation  (C-20)  then  becomes 


Equation  (C-22)  gives  the  rate  of  mean  Rothalpy  variation 
along  the  mean  meridional  streamline 

For  the  axisymmetric  case,  _  yv  ^  X 

3>t  ~ 

and  according  to  the  definitions 


(C-23) 


(024 ) 
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APPENDIX  D 


DERIVATION  OF  MEAN  FLOW  EQUATIONS  OF  MOTION  FOR 
STREAMLINE  CURVATURE  TECHNIQUE 

Equation  of  Motion  in  1-direction 

In  deriving  the  mean  flow  (pitchwise  averaged)  equation 
of  motion,  it  is  assumed  that  pitchwise  variation  of  density 
is  small.  Hence  the  mean  flow  equations  of  motion  derived  in 
Appendix  A  are  considered  valid  here,  variation  of  density 
in  the  streamwise  direction  is  taken  into  account  through 
the  continuity  equation. 

For  Streamline  Curvature  technique,  a  single  equation 
of  motion  is  derived  in  the  meridional  plan  (or  r-z  plane) . 
The  mean  flow  equation  of  motion  will  be  derived  in  a  direc¬ 
tion  inclined  at  an  angle  Y  to  the  radial  direction.  Figure 
(4.1)  gives  the  details  of  various  directions  and  angles  to 
be  used  in  the  derivation. 

The  procedure  adopted  for  getting  the  desired  equation 
is  to  start  with  mean  flow  equations  of  motion  in  r,  9,  and  z 
directions  and  derive  a  single  vector  mean  flow  equation  in 
Crocco's  form.  Equation  in  any  particular  direction  can 
then  be  derived  from  the  crocco  equation  by  taking  its 
scalar  product  with  the  unit  vector  in  the  desired  direction. 

Pitchwise  averaged  equations  of  motion  in  the  r,  9,  and 
z  directions,  in  a  coordinate  system  fixed  to  the  blades, 
were  derived  in  Appendix  A. 
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* 


These  are, 
r-direction 

Kl "i  +  -2"Ws-f 1--&1 


l5T+  Z3H 


Q-direction 


_  _  J_  ^  (D-l) 


^_VVe  +  vv,  ^Ws  +  2  J2.  W-t  +  VV©_^fc 

Ok  -»  *V 


—  2^  -t-  *t  ^~©T> 


(D-2 ) 


z-direct ion 


W-t  -+-W, _ I2fe  7  +  C  .  p 


where  2*.^  ,  and  2^^  are  t^ie  apparent  stress  terms 

Z>v  i dQ'^4')-  >«??  ) 

Ze,  _  I  C 2<?  *>9^  £ >  w»0*'f)  ,D.3a) 


and  ^  f  ^2  (2  Mt 


bwlwii 

DZ  *  > 


/vfc-  /•©£  and  are  the  r,  Q,  and  z  components  of  the 

mean  inviscid  blade  force,  given  by 

£.  =  -  AtP  --bo~U  j  F0,  -  -  ^fe_  cw  A>v 

^  ?>t^©  '  ?*A©  ^  'CAQ 

where  /•**.  is  the  mean  blade  lean  angle  and  is  the  mean 

blade  angle  with  respect  to  z  axis  in  the  azimuthal  plane. 
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r 


and  are  the  components  of  mean  viscous 

drag  force  in  r,  0,  and  z  directions. 


By  adding  and  subtracting  W0 
( D- 1 )  ^  Ohfe 


% 

9* 


to  equation 


to  equation  (D-3),  the  equations 
(D  1)#  (d  2) f  and  (d-3)  can  be  rearranged  as 


2-iC*£+1£+'*k)  i  g-r  ^  — 

-  vv|  (_  'dQ*wo)J  +  +  2  JlTvt 

c)  (1? 


a) 


b) 


(D-4) 


(D-5) 


/'sr 

“  “  ^  <D-6> 

Now  replacing^  by  ^  _  ^_Q  in  and  terms> 

and  since  ^  and  =  <1^  these  equations 


become 


b) 


c) 


_  =  _1  |E  -2*>  +  F*  +TW 

-3(-  +  5  f2£^))  y  C  ■= 

'  ?z  y  *  1 22.  y-  -  <4;y  ^©i  +  ^ 


(D-7) 


'0U  (D-8) 


~  2a* +&£+£*» 


( D—  9 ) 
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The  bracketed  terms  of  equations  (D-7) , (D-8)  and  (D-9) 
are  the  components  of  vorticity.  Thus  for  pitchwise  average 


flow 


(vxx)  =K" 


9H 


NxV)  _  (JM*  - 

V  Jq  -  Kdz  'Jt/ 
(vxY  )  -l[  '$(**>•)  1 


(D-10) 


We  can  now  get  a  single  equation  in  vector  form  replacing 
equations  (D-7) ,  (D-8)  and  (D-9) 

-  vv  *  (fv*V )  +  V  (£  <tf)  -  v 

~  ~  ^  -  ('t.Zu  -+-  QZe^Z-Zz)  +  -*-£o^  (D_n) 


sre  is  the  mean  inviscid  blade  force  vector  and  /v^J. 


is  the  mean  viscous  force  vector. 


We  know  from  Gibbs  relation  that 


-L  Yb  =  VA  -  T  VS 


(D-12) 


Taking  pitchwise  average,  this  becomes 

J-  Vjb  -  V£  -  T  YJ 


(D-13) 
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be 


Here  the  temperature  fluctuations  have  been  assumed  to 
small 

Equation  (D-ll)  then  becomes. 


w  x(Vx^)=  V -TVS 

+(*t2\l  +  9  Zqx  Ey  (D-14) 


Prom  the  definition  of  Rothalpy  (relative  total  enthalpy)  i, 
we  know 


X  ~  A  -h  SI.  /Z$g 

2. 


(D-15) 


Using  the  definition  (1.1) 


for  X 


we  have 


_  2 


I  =  k  -t  ~f  /z  -'tJIVv 

—  __2  _ 

Substituting  for  (  fo  -f-  — ‘t-Q.'Tfa  )  in  (D-14),  we  get 

vv  x  (y  xX)  =  VJ  -  T  vJ 

-f  &Z&>  -fty-Rf  j  (D-16 ) 


Equation  (D-16)  is  the  pitchwise  averaged  counterpart  of  the 
steady  state  Crocco  equation  modified  to  include  viscous 
terms. 


This  is  the  basic  equation  for  deriving  mean  flow  equa¬ 
tion  of  motion  in  any  desired  direction. 

in  order  to  arrive  at  the  mean  flow  equation  of  motion 
in  1-direction  (Figure  4.1),  we  will  take  dot  product  of 

A 

equation  (D-16)  with  unit  vector  1. 
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Let  /A  =  VV  X  C  V  X  V^_  ),  and  let  At ,  >Ae  cltlcL  Ag 
be  the  components  of  vector  A  in  r,  9  and  2  directions. 

Let  j  F-ni  and  Fq^  be  the  components  of  invis- 

cid  force  A'  wt^’iand  O  directions,  and  Fiv  and 

F&v  be  the  components  of  viscous  force  in  m,  n  and  9 

/A 

directions.  Writing  the  equation  for  dot  product  /  • 

l  •  It  A*  +eA0+zA£1__  'zT-  t2L 
l  J "  7T  * 


*  r  A  A  A  ") 

(D-17) 

-(•{**  + 


Since  /  lies  in  the  meridional  plane  which  is  per- 

A 

pendicular  to  9  direction  L^Q-O 

From  Figure  (4.1)  we  find  .-c'^fc  -  C&J  T 

l.  z  =  «/cn, T  t  / '  m  =  x£tr^  (yf  -&r)  j  /»-??-  Cod(F+r) 

Here  T  is  fixed  and  ^5  is  defined  as  the  angle  of  mean 
stream  surface. 

Substituting  these  relations  in  equation  (D-17),  we  get 

£  v£ a  (vx X )J^ O r  +  ^ x (y xV)J^ r 
=  It  -_r|/  +  ZuCvr* 

W  *1  _  (D-18 

—  i  +  FyKl ‘dvn  c& ~  ( Eni  +Fnv)£to(?'hr) 
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Simultaneous  solution  of  these  two  equations  gives 


3D  ) 

—  _ oTtl _ _ 

Jusri,  <$ALix  T  —  sCoJ  pSsC&d  T 


(D-20) 


2()_  Jonl^iL  sC*r£S^ 


(D-21) 


We  can  now  transform  Terms  i  and  n  in  the  left  hand 
side  of  equation  (D-19)  in  terms  of  derivatives  in  m  and  1 
directions. 

Using  equations  (D-20)  and  (D-21),  Term  I  becomes. 


Term  I 


=  f  [-  +  CWfc*  #tWej 

^  r-XoAf! Ccur  Jim.jt/A».r-£*4pG>dr 


Simplifying 


Term  I  — 


_  hg  f 


*  u 


Term 


changing  to  W©  +  X.SL  ,  We  get 


Substituting  in  (D-18)  for  andjy£*fox^)J 

om  the  following  relations 

|  ]V*  (yx  J)j  _  fa  ^  V&  )  ,  Wa  Q  ftU.; 


(Note : 

and  rearranging,  gives 


W®  (■CojrJ^KJinr^^f 

^  ( _ _ J 

(Zl  Z  " 

+  ‘  (WiAunT  -  'IZz&J'f)  J  =  21  _T3S 


(D-19) 


■b  (Gur)2jc^  +(di n  r)2Z),  -  ( frn.j4'  ^v)Acn  (f+r)  - 

In  order  to  get  radial  and  axial  derivative  in  terms  of 
derivatives  in  m  and  1  directions,  following  transformation 
relations  are  used 


2_o 

Irry^ 

_  2fQz 

*  .20  2i 

9't  9»t 

3D 

DU  9t 

sr 

o|? 

3D 

9?*t 

=  4od^^J*f 

Aonfl  2S ) 
'b\. 

90 

-  JinT^U 

+  'Cur— 

(D-19a) 
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or 


Term  I  -  ^9  2^We)  +  2  W6  JL  Gj 

x  at 


(D-22 


For  evaluating  Term  II ,  we  use  the  relations 
~  'don ft  Cv*vd,  'Vb  =  ^nGd  $ 


This  gives 


Term  n 


=  (^iuri /  din  T  -  CdA  rj 


Using  equations  (D-20)  and  (D-21),  we  get 


ZL  - 


dvnr-  C&i /  Gjr)  X 
-  CxrWt-J-  (dcnr^ 


-  ^^5, 


=  %>  (Ai^JonT _  CojyiCajrJ  x 
I  ^don*  0  - G<j/  Zad^ 


f 


_ a  _2  _ 

Simplifying  and  using  the  relations  -^h 

_/  ,  we  get 

9m” 


_  _ a 

■Tern  II  =  ^  Gj  (jfi+rj  (D_23) 

11  ***>  \ 


Substituting  Term  I  and  Term  II  in  equation  (D-19) , 

rearranging  the  equation  and  using  the  relation  _1_ 

Gjy  -f-  'Aur\T~  •  equation  (D-19)  becomes 

“3*.  1Z 


Aon(0+r)  -h  tyn  'CyiCft+r)  _  2  fa  Wo) 

-  ^  -  2Wq  SLCzjT  _  t*2J 
^  3/ 


'Cojrf±  ift *n(i)  ?(>wi*gr)  w'?/ 
v-  /  t  31  y 


•f 


>f  /  a*  -az  +  —  / 


knr 
/f  /  : 

—  ,)  (0+r) 
-  ( Fn  i  ■+■  Frnf)  CaAQt+r) 


(D-24 ) 
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This  is  the  mean  flow  equation  for  Streamline  Curvature  Computa¬ 
tional  method.  This  equation  has  to  be  suitably  modified  for 
the  computer  program  depending  on  the  type  of  input  data. 
Inviscid  and  viscous  force  terms  have  also  to  be  expressed  in 
terms  of  the  known  input  data. 

Mean  Flow  Equation  of  Motion  in  the  Meridional  Direction 

Pitchwise  averaged  equations  of  motion  (D-l)  and  ( D— 3 )  in 
the  radial  and  axial  direction  can  be  written  as 


9*  £  ?Z.  X  J  ^ 


(D-25) 


b)  W*.  3WZ 


( D-26 ) 


0  i  —  _  ^ 

Here  (  «+■  <^,)  has  been  replaced 

and  =■  +  Fvo  j  ^z.  ~  F&- 

X 

Using  the  relations 

"S;2=  =  Wz  +  Wt  ; 

=  "Sin  $  j 

and  after  some  algebraic  manipulation,  we  get 


dm  -t 


= - ^  r}_£  +  +  f^G><5<f>) -  (D-27 ; 


Frn  -  f-'L  dvrirf  -f- 


«MWM 


This  reduces  equation  (D-27)  to  the  following  form 

T  Hn  =  ^  -  (TtxJL/  *2^  <^) 


'D-28) 


is  eliminated  from  equation  (D-28)  by  using 


equation  (40) 


4  = 


-t  -a.* 


Zast<x».ji  - t  2 ^ 


Equation  (D-28)  becomes  for  the  bladed  region 

-  6x1/)  _  2^ 


(D-29) 


Equation  (D-29)  is  the  mean  flow  momentum  equation  in  the 


meridional  direction. 

Now  using  the  equation  ~T~cLS -  (j?dT  -  d-J* 
From  the  equation  of  state  7=  ,  we  have 


(D-30) 


dT-z(f  -ft) 


Using  this  along  with  the  relation  E  -  Zd  ,  in  eoua- 

Q>  "  r 

tion  (D-30) ,  and  taking  derivative  with  respect  to  m,  we  get 

fi f- 

-  a  —  *■ 

Taking  Zt  -Q,  - 

-  f  7C 

and  substituting  this  for  y  "C*  ,  in  equation  (D-29),  we  get 

r  a<>H. 

the  modified  form  of  meridional  direction  equation  of  motion: 
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Figure  1.  Rotor  section  cut  at  z  =  constant  plane. 
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Figure  3.  Meridional  plane  picture  of  a  passage. 


Figure  S,  Apparent  shear  stress  v'v'  ,  0.1  chord  downstream 

OF  THE  ROTOR. 


25000 


Pressure  Ratio 


Figure  21.  pg  /pt  ,  0.1  and  1.0  chord  downstream  of  the  rotor, 
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START 

I 

READ  INPUT  DATA 

ESTIMATE  STREAMLINE  PATTERN 

j!PASS=H 


c^n 

INTERPOLATE  INLET  CONDITIONS 
IF(IPASS.GT.l)  » 


ESTIMATE  'INLET  VELOCITY 

SET  VMOLD( J)=VM( J , I) 


|SET  ENTROPY , ENTHALPY  ETC.' 
X 


IF(IPASS.CT.l)  "  ■» 
SET  VM(J,I)=VM(J,  1-1)1 


VMA 


CF (STATION  LIES  OUTSIDE  THE  BLADED  REGION 

INTERPOLATE  STATION  DATA 
DETERMINE  ENTROPY,  ROTHALPY  OR 


ENTHALPY  AND  WHIRL  ANGLE 


SET  VMLOLD(J)=VM(J, I) 
flTER=l| 

\ 


SOLVE  MOMENTUM  EQUATION 
ESTIMATE  MASS  FLOW  RATE  W 
COMPARE  W  WITH  SPECIFIED  MASS  FLOW 

IF«j.K.)* 


RESET  V  (MIDSTREAM) 

I  m 

TER=ITER+11 


COMPARE  V  WITH  V  (IPASS-1) 

•  m  m 

RE-ESTIMATE  STREAMLINE  LOCATIONS  AND  STORE  THE 
IF(I=LAST)— ,  CHANGES  REQUIRED 


r  i 

I  =  1+1 


IF(VM.EQ.VM(IPASS=1))  FOR  ALL 

- ONS" 


SET  fcpASS=IPASsSl 


RELOCATE  THE  STREAMLINES  BASED 
ON  THE  STORED  CHANGES 

r~-  ■  - 


STOP 


IPASS  PASS  NUMBER 

I  STATION  NUMBER 

V  (J,I)  MERIDIONAL  VELOCITY  FOR  THE  CURRENT  PASS 

AT  THE  Ith  STATION  AND  Jth  STREAMLINE 
ITER  ITERATION  NUMBER  FOR  ITERATIVE  SOLUTION  OF 

MOMENTUM  AND  CONTINUITY  EQUATIONS 

Figure  2‘\.  Overall  logic  of  the  original  streamline-curvature  program. 
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OVERLAY  UD0300  J  SUBROUTINE  CORLTA  SUBROUTINE  SPLINE 

IF(NCOR.EQ.l)  CALL  SUB-CORLTA  j w Reads  Lata  on  Velo-  .  Calculates  Cubic 

'  city  Correlations  Splinefit  Coeff- 

I  icients 


main  program 
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1)  Mass  Plow  Rate  =  81.88  lb  /sec 

m 

2)  R.P.M.  =  13193 

3)  Total  Pressure  (P^)  at  inlet  =  14.7  psia 

4)  Total  Temperature  (T^)  at  inlet  =  520.0°  R 


Relative  Flow  Angles 


Station  4 

Station  5  ! 

Station  6 

Radial 

Location 

(inches) 

Relative 

Plow 

Angles 

(degrees) 

Radial 

Location 

(inches) 

Relative 

Plow 

Angles 

(degrees) 

Radial 

Location 

(inches) 

Relative 

Plow 

Angles 

(degrees) 

7.4560 

-51.700 

7.4560 

-35.576 

7.456 

-23.000 

7.9818 

-54.720 

7.9820 

-43.441 

8.2977 

-36.250 

9.4365 

-59.800 

9.4365 

-53.500 

9.3800 

-49.250 

10.577 

-62.240 

10.577 

-58.720 

10.140 

-52.700 

11.650 

-65.140 

11.650 

-65.200 

11.650 

-65.000 

Figure  28a.  Basic  input  data  for  computation. 
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Station  6 
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Units  for  Velocity  Correlations  Ft2-Sec"‘ 


Figure  28b,  Apparent  stress  data, 
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Experimental 
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r/rt 

Figure  37.  Apparent  stresses  (Vg/r),  0.1  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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'  — nnrflWMnii 


Experimental 


Calculated  with 
Cascade  Losses  and 
Apparent  Entropy 
(Station  6) 


_l _ 1 _ I _ 

°«7  R/Rt  °*8  °«9 

Figure  Ml,  Mean  pitchwise  velocity,  0,1  chord 

DOWNSTREAM  OF  THE  ROTOR. 


0.5  0.6  0.7  R/R^  0.8  0.9  1.0 

Figure  42,  Mean  axial  and  radial  velocities, 

0,1  CHORD  DOWNSTREAM  OF  THE  ROTOR. 
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Figure  47.  Mean  axial  and  radial  velocities, 

0.1  CHORD  DOWNSTREAM  OF  THE  ROTOR. 


187 


V2  Ft  Sec 


N? 


Experimental 


'Calculated  with 
Cascade  Losses  and 
Rothalpy  variation 
(Station  7) 


Figure  49.  Mean  axial  velocity,  1.0  chord 

DOWNSTREAM  OF  THE  ROTOR. 


Figure  50.  Mean  radial  velocity,  1,0  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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Figure  52.  Mean  pitchwise  velocity,  0.1  chord 

DOWNSTREAM  OF  THE  ROTOR, 
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Figure  54.  Mean  pitchwise  velocity,  1.0  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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Figure  57.  Mean  pitchwise  velocity,  0.1  chord 

DOWNSTREAM  OF  THE  ROTOR. 


Figure  58,  Mean  axial  and  radial  velocities, 

0.1  CHORD  DOWNSTREAM  OF  THE  ROTOR. 
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Mean  axial  velocity,  1.0  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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Rothalpy  Variation 
(Station  6) 
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Figure  62.  Mean  pitchwise  velocity,  0.1  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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Figure  64.  Mean  axial  velocity,  0.1  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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Figure  66.  Mean  pitchwise  velocity,  1.0  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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Figure  68,  Mean  radial  velocity,  1,0  chord 

DOWNSTREAM  OF  THE  ROTOR. 
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